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^ ABSTRACT 
> 

Context. Our ability to extract information from the spectra of stars depends on reliable models of stellar atmospheres and appropriate 
techniques for spectral synthesis. Various model codes and strategies for the analysis of stellar spectra are available today. 
Aims. We aim to compare the results of deriving stellar parameters using different atmosphere models and different analysis strategies. 
The focus is set on high-resolution spectroscopy of cool giant stars. 

Methods. Spectra representing four cool giant stars were made available to various groups and individuals working in the area of 
spectral synthesis, asking them to derive stellar parameters from the data provided. The results were discussed at a workshop in 
Vienna in 2010. Most of the major codes currently used in the astronomical community for analyses of stellar spectra were included 
in this experiment. 

• • Results. We present the results from the different groups, as well as an additional experiment comparing the synthetic spectra produced 

^ by various codes for a given set of stellar parameters. Similarities and differences of the results are discussed. 

Conclusions. Several valid approaches to analyze a given spectrum of a star result in quite a wide range of solutions. The main causes 
for the differences in parameters derived by different groups seem to lie in the physical input data and in the details of the analysis 
method. This clearly shows how far from a definitive abundance analysis we still are. 

Key words. Stars: atmospheres - Stars: late-type - Methods: analytical - Stars: fundamental parameters 
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* Based on observations obtained at the Bernard Lyot Telescope 
(TBL, Pic du Midi, France) of the Midi-Pyrenees Observatory, which 
is operated by the Institut National des Sciences de I'Univers of the 
Centre National de la Recherche Scientifique of France. 
** Tables 6 to 11 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-H A/| 



1. Introduction 

Spectroscopy is the basic tool of modem astrophysics. It is the 
key for revealing the elemental composition and the physical 
conditions in the spectrum-forming layers of stars. Interpreting 
the information contained in the spectra requires knowledge 
about the physics of the stellar atmosphere, the line formation 
processes, and the atomic and molecular data. Parameters de- 
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rived from the analysis of high-resolution spectra via compari- 
son with stellar models have great potential but suffer from sys- 
tematic uncertainties due to insufficient input physics in both the 
model atmospheres and the spectral synthesis, as well as to dif- 
ferent fitting approaches. 

Several codes to calculate atmospheric models exist today 
and are used by various groups around the world to analyze both 
spectroscopic and photometric data. However, the implementa- 
tion of the physics, the atomic and molecular data used, and the 
details of the method of deriving stellar parameters from the ob- 
served data differ among the various research groups. Therefore, 
a comparison of the various codes and their output may help us 
understand uncertainties in the analysis of stellar spectra intro- 
duced by these various components involved in the fitting pro- 
cess. These uncertainties also have a major impact on the inter- 
pretation of photometric data or the modelling of stellar popula- 
tions. 

In this paper we present a comparison of a variety of model 
codes that attempt to analyze the spectra of cool giants. Red gi- 
ant stars are quite challenging targets for modelling with their 
complex atmospheres and the large number of lines, in particular 
those of molecular origin. Hence, they provide a good testbed for 
exploring the validity of input physics, line data, and modelling 
approaches. The aim was to test how comparable the results are 
when applying different methods. 

For this comparison, three of us (U. Heiter, T. Lebzelter, W. 
Nowotny) designed the following experiment, hereafter referred 
to as Experiment 1: colleagues, who are engaged in stellar pa- 
rameter and abundance determinations for cool giants on a regu- 
lar basis, received high resolution and high signal-to-noise ratio 
(S/N) spectra of four cool stars and were asked to derive their 
fundamental stellar parameters, such as the effective temperature 
(Teff), the surface gravity (log g), and the metallicity ([Fe/H 
In addition, we provided corresponding photometric data in vari- 
ous bands. However, no identifications of the sources were given 
in order to prevent the participants from comparing their findings 
with data in the literature. The list of authors of this paper illus- 
trates the high level of interest in this experiment. The results 
were compared and discussed during a workshof|^ held at the 
University of Vienna, August 23-24, 2010. 

The choice of cool giants as targets was driven by the moti- 
vation of this experiment within the framework of ES A's upcom- 
ing Gaia missiorj^ In preparation for the exploitation of a large 
amount of spectroscopic and photometric data with the aim of 
determining accurate stellar parameters, there is a clear need to 
identify key areas where model spectra can and should be im- 
proved, and to determine the influence of different methods of 
analysis on the results. Giant stars will play an important role 
within the sample of objects that will be studied by Gaia. 

At the workshop we agreed to perform a second comparison 
of our models, hereafter referred to as Experiment 2. In this case, 
each participating group was asked to calculate a high-resolution 
model spectrum in a pre-defined wavelength range using a given 
set of stellar parameters. 

This article is organized as follows. In Section [2] we pro- 
vide details on the design of the experiments, and summarize the 
properties of the two benchmark stars analysed in Experiment 1, 
a Tau and a Cet. A detailed description of the modelling ap- 
proaches used in the experiments is presented in Section [3] In 



Sections]?] and |5] the results of both experiments are presented 
and discussed. As a by-product we give revised stellar param- 
eters for the two benchmark stars. The outcome of our experi- 
ments forms a basis for future improvements of stellar spectrum 
modelling. This will be an important step towards accurate stel- 
lar parameters of giant stars observed by Gaia and Gaia follow- 
up programmes. 



2. Experiment set-up and benchmark stars 

2.1. Experiment 1 - stellar parameter determination 

High-resolution and high S/N spectra of four targets were pro- 
vided for Experiment 1. Two of them were of real stars, namely 
a Tau and a Cet, and covered the visual range of the spectrum 
between 4900 and 9750 A. They were obtained by one of us (U. 
Heiter) using the NARVAL spectrograph at the 2m Telescope 
Bernard Lyot atop Pic du Midi (Auriere 2003). The resolution 
was set to 7;=80000. S/N>200 was achieved throughout the 
whole spectral range for both stars. The data were reduced with 
the Libre-ESpRIT pipeline ponati et aL]|1997| ). The extracted 
and calibrated echelle orders were merged by cutting the orders 
at the centers of the overlap regions. The spectra were not cor- 
rected for telluric features, but a spectrum of a telluric standard 
star taken in the same night was provided. Both objects have 
been studied in detail in the past and have been used as refer- 
ence targets in several investigations. In Table [T] we summarize 
the stellar parameters which we assume are the most accurate 
ones available in the literature. A more detailed description of 
the targets is given in Sections |2.2| and |2.3| Three wavelength 
ranges were recommended, on which the experiment should fo- 
cus: [4900 - 5400], [6100 - 6800] and [8400 - 8900] A. Only 
the second of these regions was judged to be significantly con- 
taminated by telluric lines. 

The other two spectra were synthetic ones computed for re- 
alistic stellar parameters. Within the experiment they should al- 
low for a more direct comparison between the models with- 
out the uncertainties of stellar parameters and the unidentified 
features we see in observed data. For the artificial data we 
used two COMARCS model spectra (for a description of the 
COMARCS models see Sect.|3.1.1[) calculated by W. Nowotny 



' [Fe/H]= log ^ - log [M/H] is defined accordingly using 

any atoms heavier than He. 
^ Kindly funded by the ESF within the GREAT network initiative. 
^ http://sci.esa.int/gaia/ 



and T. Lebzelter. To simulate observational effects some gaus- 
sian noise was added (S/N=125) and the output was rebinned 
from 7?=300 000 to R-50 000. The chosen model parameters are 
listed in Table [2| The input parameters for Star 3 were chosen to 
reproduce a slightly metal poor K-type giant as might be found 
in an Large Magellanic Cloud cluster. Star 4 should resemble a 
typical field star on the Asymptotic Giant Branch (AGB). The 
corresponding quanti ties follow the predic tions from the stellar 
evolution models by Marigo et al. (2008 1. The ratio of '^C to 
'^C was solar in both cases. For the synthetic spectra we chose 
a wavelength range of A =15 456-15 674 A. Within this part of 
the spectrum one finds lines of CO, OH and CN as well as sev- 
eral atomic lines. The wavelength range is almost free of tel- 
luric lines (e.g. Hinkle et al.||1995) . The synthetic spectra have 
a much smaller wavelength coverage than the NARVAL spec- 
tra which reflects the fact that today's near infrared spectro- 
graphs that reach a resolution of i?=50000 (CRIRES, Phoenix) 
also cover only a comparably small wavelength range at a time. 
Broad band colours of the model stars, listed in Table]2]as well, 
have been calculated from lower resolution spectra (/?=10000) 
over a wavelength range from 0.45 to 2.6 yum as described in 
[Nowotny et al.| ( |2i31 l| l. The participants of the experiment were 
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Table 1. Stellar parameters of a Tau (Section 2.2 1 and a Cet 



(Section 2.3 1 based on published data. 
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Ref. 


a Cet 


Ref. 
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Star 2 




Name 




Aldebaran 




Menkar 




riK 




14j / 




oil 




MK type 




K5 III 




Ml. 5 nia 




^ eit \^^/ 




3930±40 


fa) 


3800+60 


(al 

V"/ 


T,ff (K) 

^ ett V^*-/ 




3920±130 




3730±75 


fi) 


L [Lq] 




440±20 


(c) 


1870±130 


(c) 


log (g (cm 




1.2±0.1 


(a) 


0.9±0.1 


(a) 


log {g (cms" 




1.2±0.5 


(b) 


0.7±0.3 


(i) 


M[Mo] 




1.3±0.3 


(d) 


3.0±0.5 


(d) 


[Fe/H] 




-0.22±0.11 


(e) 


+0.02±0.03 


(i) 


V-/ 




2.17±0.02 


(f) 


2.51±0.02 


(f) 


J-K 




0.97±0.03 


(f) 


1.08±0.03 


(f) 


V-K 




3.67±0.03 


(f) 


4.21±0.03 


(f) 


fmd [km s-'] 




54.26±0.03 


(g) 


-26.08±0.02 


(g) 


J) sin I [km s 




5±1 


(h) 


3±2 


Q) 



Notes. direct parameters (see text). mean of parameters used for 
spectroscopic [Fe/H] (see Tablejsji. from bolometric flux and parallax 
(see text). from T^.^ and L, and evolutionary tracks (see text). see 
Table [3] Johnson et al. (1966 |. heliocentric ra di al velocity from 



Melendez et al. 



Famaey et al.,l,2005j. "",Hekker & Melendez| l |2007^ 
p.QOS\ . w |ZMnanov et aLn2008> . 



Table 2. Input parameters and model colours for the two artifi- 
cial stars in Experiment 1 . 





Star 3 


Star 4 


Teff (K) 


4257 


3280 




319 


3816 


log (g (cm s"^)) 


1.47 


0.06 


M[Mo] 


1.165 


1.509 


[Fe/H] 


-0.4 


+0.1 


C/O 


0.35 


0.55 


V-I 


1.25 


3.58 


J-K 


0.82 


1.23 


V-K 


2.94 


6.89 



not informed about the artificiality of these data. All four spectra 
used in Experiment 1 are available onlin^ 

A specific list of atomic line data was provided, and was 
suggested, but not required, to be used for the analysis. This 
line list was extracted on 2010-04-29 fro m the Vienna Atomic 
Line Database (VALD) (Uppsala mirro:|^ Piskunov et al. 1 1995 



Kupka et al.|1999| [Heiter et al. 2008 ), and included all lines in 



the database within the three recommended optical wavelength 
regions and the single infrared region. 

2.2. a Tau 

Stellar surface parameters (Tetj and log g) can be determined ei- 
ther from angular diameter measurements in combination with 
additional data (called direct parameters hereafter), or from a 
model atmosphere analysis of photometric or spectroscopic data. 
The direct Teff value is obtained from the angular diameter 6 
and the bolometric flux Fboi according to Eq. [T] where cr is 
the Stefan-Boltzmann constant. The direct log g value is derived 

ftp://ftp.astro.uu.se/pub/Spectra/ulrike/ComparativeModelling/ 
' http://www.astro.uu.se/<tilde>vald/php/vald.php 



from 9, the stellar mass M, and the parallax n according to Eq.|2] 
where R is the linear stellar radius and G is the constant of grav- 
ity. 



Fhoi 



' eff 



o-i0.59f 



GM 

g = — R 



In 



(1) 
(2) 



The angular diameter of a Tau was determined recently by 
|Richichi & Roccatagliata| ( |2005[ ), using both lunar occultations 
and long-baseline interferometry (VLTI- VINCI, A'-band), and 
taking into account limb darkening. The integrated absolute flux 
was measured for a Tau b y |di Benedetto & Rabbia| ( |1987] ) and 
|]VIozurkewich et al.| ( |2003] l. The measured Q = 20.58 + 0.03 mas 
and Fboi = 33.57 + 1.35 nW m"^ result in the direct T^s value 
given in Table [T] With the H ipparcos parallax n - 48.92 ± 
0.77 mas ( |van Leeuwen|2007| l, this results in 7? = 45.2 + 0.7 
and in the luminosity L given in Table [T] We estimate the mass 
of a Tau using two different sets of stellar evolutionary tracks, 
those pubUshed by the Padova group (Bertelh et al.l2008 2009^ 
and the Yonsei-Yale (Y2) models ( |Yi et al. 2003; Demarqiie] 



|et al. |2004| l. For solar metallicity tracks (Padova: Z=0.017, 
F=0.26; Y2: Z=0.02, X=0.71), the direct T^^ and L imply a 
mass of 1 .6/V/0 for both model grids, while the metal-poor tracks 
corresponding to [Fe/H]=-0.3 (Padova: Z=0.008, 7=0.26; Y2: 
Z=0.01, X=0.74) suggest a mass of 1.0 to 1.1 M^. Interpolating 
between these values results in the mass in Table [T] Thus, we 
arrive at the log g value for a Tau given in Table [T] 

a Tau has been studied with high resolution, high S/N spec- 
tra in nine publications since 1980 (according to the PASTEL 
catalogue, |Soubiran et al.|2010| l. The stellar parameters derived 
in these works and the references are given in Table [3] Effective 
temperatures used in the spectroscopic analyses have been de- 
rived from various photometric calibrations by most authors. 
Combining the results given in six publications, the mean pub- 
lished photometric Teft of a Tau is 3850±40 K. This is in good 
agreement with the latest value of Teif =3880±40 K determined 



with the infrared flux method (IRFM), by Ramirez & Melendez 
( |2005 1, which is an update of the work by |AIonso et al. \ 1999| l 



Two publications from 1981 and 2007 derive Teff in a spectro- 
scopic way (excitation equilibrium of iron line abundances) from 
high-resolution spectra in the optical wavelength range. They ar- 
rive at almost the same value, close to 4120 K, which is sig- 
nificantly higher than the photometric one. On the other hand. 



Melendez et al. (2008 1 obtained a spectroscopic red of 3890 K, 
close to the IRFM value, based on high-resolution infrared spec- 
tra centered on 15555 A. 

The method for determining the surface gravity varies sig- 
nificantly between the publications. Four of the authors derive 
loggf in a spectroscopic way in the optical wavelength range 
(ionization equilibrium of iron line abundances) and arrive at 
a mean value of 1.1+0.5. Two authors use absolute magnitudes 
and stellar evolution calculations, and cite a higher mean value 
of 1.6±0.1. The spectroscopic log 3 of Melendez et al. ( |2008| l 
determined from IR spectra is close to the values obtained from 
optical spectra. The highest log g value of 1.8 is determined from 
DDO photometry by [Fernandez- Villacanas et alT| (jT990}), who 
cite an error of 0.2. 

The metallicity of a Tau is determined in eight publications 
and found to be below solar (mean value -0.2+0.1, see Table[T]i. 
The results can be divided into two groups, three authors us- 
ing reff>3900 K and five authors using reff<3900 K. The cor- 
responding [Fe/H] values cluster within a few tenths of a dex 
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Table 3. Results and references for previous spectroscopic 
works on a Tau. 

Teff (a) log 9 (b) [Fe/H] Reference 



4140 


sp 


1.0 


sp 


-0.33 


3830 


ph 


1.2 


sp 


-0.14 


3850 


ph 


1.5 


ev 


(c) 


3800 


ph 


1.8 


ph 


-0.17 


3910 


ph 


1.6 


ev 


-0.34 


3875 


ph 


0.6 


sp 


-0.16 


3850 


ph 


0.6 


li 


-0.10 


4100 


sp 


1.7 


sp 


-0.36 


3890 


sp'' 


1.2 


sp'' 


-0.15 



Lambert & Ries^ (|1981j 
Kovacs ( ■ " 
Smith & 



ert 



1 (11985 I 



hernandez-ViUacanas et al. 
McWilliam(1990, 



(19901 



Luck & Challenerl 
Mallik( 19981 



19951 



Hekker & Melendez d2007 1 



Melendez et al. i 200 



Notes. Method for T^ff determination: sp . . . spectroscopic (excita- 
tion equilibrium of iron line abundances), ph . . . photometric calibra- 
tions. **' Method for log g determination: sp . . . spectroscopic (ioniza- 
tion equilibrium of iron line abundances), ph . . . photometric calibra- 
tions, ev . . . from absolute magnitudes and stellar evolution calcula- 
tions, li . . . from literature. a Tau was used as a reference object for 
the spectroscopic analysis of M giants. *''' from infrared spectra. 



around -0.35 and -0.15, respectively. The mean Teff and \ogg 
values used for deriving the mean [Fe/H] are given in Table [T] 

The spectral type K5 III given for a Tau in Table [T| was first 
published in the original MKK Atlas ( [Morgan et al.| 19431 ) and 
has been quoted throughout the literature. However, the star's 
TiO strength, as measured by narrow-band classification pho- 
tometry ( )Wing||2011) , yields a spectral type of K5.7ni on a 
scale where type MO.O immediately follows K5.9. The photo- 
metric colours also indicate an effective temperature closer to 
MO.O than to K5.0. 



2.3. a Cet 



The angular diameter of a Cet was determined by Wittkowski 



et al. (2006a I from long-baseline interferometry, using the same 
instrumentation as for a Tau, and taking into account limb dark- 
ening. The same authors also determined Fboi from integrated 
absolute flux measurements. The measured 6 = 12.20+0.04 mas 
and Fboi = 10.3 ± 0.7 nW m"^ result in the direct T^ff value given 
in Table [T] With the Hipparcos parallax n = 13.10 ± 0.44 mas 
( [van Leeuwen|2007|l, this results in R = 100 + 3 and in the 
luminosity L given in Table [T] Using this luminosity and the di- 
rect Teff value, we estimate the mass of a Cet (see Table [T i from 
evolutionary tracks for solar metallicity (see Section 2.2 1. The 
masses from the two different sets of models agree within 0. 1 Mq 
The 2007 Hipparcos parallax is smaller than the "original" one 
by about 10%, which results in a mass 30% higher than derived 
by |Wittkowski et aL]p006a| ). The mass estimate could be refined 
in the future through asteroseimology, but in any case the uncer- 
tainties for the direct log g value are already much smaller than 
the spectroscopic uncertainties. As for a Tau, from these stellar 
data we can derive an accurate direct log g value for a Cet, given 
in Table [U 

For a Cet, there are no previously published high-resolution 
spectroscopic studies in the optical wavelength range. The T^g 
value determined with the IRFM value is 3720+50 K ( ,Ramfrez| 
& Melendez 2005 |. The star is in cluded in the infrared spectro- 
scopic study of Melendez et al. ( 2008), who determine a spec- 
troscopic Ted close to the IRFM value, a spectroscopic log g- 
0.7+0.3, and solar metallicity (see Table[T]i. The TiO-based spec- 
tral type of a Cet, from narrow-band classification photome- 



try ( [Wing |20TT] l, is Ml. 7 III, in substantial agreement with the 
Morgan-Keenan (MK) type shown in Table [T] 

2.4. Experiment 2 - comparison of syntlietic spectra for fixed 
parameters 

For this experiment the set of stellar parameters was prede- 
fined in order to be able to compare the output of the various 
combinations of models and spectral synthesis codes directly. 
Participants were asked to compute spectra at high resolution 
(R > 300000) for T^g =3900 K, log 3= 1.3, and [Fe/H]= -0.2, 
i.e. a parameter set close to the values corresponding to a Tau. 
Microturbulence was fixed at 2.0km s"', and a mass of 2 Mq was 
to be assumed. The synthetic spectra covered the three optical 
wavelength regions that were recommended for the analysis of 
Star 1 and Star 2 in Experiment 1. The same abundance pattern 
( Asplund et al.|2009| was used by all participating groups. 



3. Modelling 

In this section we summarize the main characteristics of the 
models used for the experiments. In our sample we have two 
major model 'families' - MARCS (Model Atmospheres in a 
Radiative Convective Scheme) and ATLAS - where several im- 
plementations and individual further developments of the origi- 
nal code were employed, and three alternative models. For the 
discussion below we introduce abbreviations for each model 
and implementation, e.g. 'M' for MARCS-based modelling, 
and 'Ml' for a specific participating team applying the code. 
Where possible we refer to more extensive descriptions of 
the codes published elsewhere. Comparisons of different atmo- 
spheric models can be found in the original literature describing 
the codes. For example, [Gustafsson et al. ( 2008[ l write on the 
comparison between MARCS and ATLAS models: 

"In view of the fact that these two grids of models are 
made with two totally independent numerical methods 
and computer codes, with independent choices of ba- 
sic data (although Kurucz's extensive lists of atomic line 
transitions are key data underlying both grids), this over- 
all agreement is both satisfactory and gratifying." 

Each group worked with different subsections of the avail- 
able wavelength range. For the optical spectra. Fig. T] shows the 
wavelength intervals used by the eleven groups who analyzed 
these spectra. Most groups analysing the spectra of Star 1 and 
Star 2 determined radial velocities, which are given in Tables 12 
and 13 and corrected the observed spectra to laboratory wave- 
lengths. The groups used different approaches to adjust the con- 
tinuum scale of these spectra, which are described in the indi- 
vidual subsections. 

The groups employing MARCS model atmospheres in gen- 
eral used the solar chemical composition from Grevesse et al. 
( |2007[ l as a reference (which was used for the calculation 
of the on-line database of MARCS models). However, Ml 



(Section 3. 1 . 1 1 adopted the abundances from Grevesse & Sauval 
(1994) for C, N, and O, and from'Anders & Gr evesse |( |1989| l for 
all other elements. A3 used only the latter source for their refer- 
ence abundance pattern. A4 (Section 3.2.4 1 took their reference 
abundance pattern from|Grevesse & Sauval ( 1998[l. The analysis 
of A5 (Section [3.1.5| l refers to the abundances of A splund et aT] 
(2005). M5, Al, and A2 ( Sections [3T5] [3XT] a nd |3.2.2| ) used 
the solar abundances from Asplund et al. ( 2009[l, and P and C 
(Sections [33] and [14| those of [Grevesse et al.| ( |1996| ). T used 
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500 550 600 650 700 750 
Wavelength [nm] 



800 



850 



900 



Fig. 1. Wavelength intervals used by the eleven different groups 



analysing the optical spectra. See §§ 3.1 to 3.5 For the M3b 



analysis (§ 3.1.31, discontinuous wavelength regi ons totaling 



"435A were used. The M6, Al, and A4 groups (§§ 3.1.6 



3.2.1 



and 3.2.4 1 used spectrum synthesis for individual spectral lines 
distributed over the indicated regions. The A2 and A5 groups 
(§§ 3.2.2 and 3.2.5 1 used equivalent widths for individual spec- 



tral lines. 



abundance values from several sources summarized in Table 4, 



case (b), of Tsuji ( 2008 1, and in Tsuji ( 2002 1. 



Some of the groups determined stellar parameters from the 
photometric colours provided together with the spectra for each 
star. These parameters, and the corresponding photometric cali- 
brations, are summarized in Tables[T2]to[T5| Table|4]summarizes 
the range of parameters of the model grids used by several of the 
groups for the determination of the best-fit spectrum. 



Sections 3.1 to |3.5| contain detailed accounts of the analysis 
work of the participating groups. In these sections. Star 1 and 
Star 2 refer to a Tau and a Cet, respectively. Each modelling de- 
scription also includes a brief discussion of the individual fitting 
results. The impatient reader may at this point skip to Section]?] 
where we start with an overview of the main aspects of each 
analysis, and the detailed descriptions may serve for later refer- 
ence. A general discussion comparing the results from the vari- 
ous groups is provided in Section ]5] 

3.1. MARCS model atmospheres (M) 
3.1.1. Ml 

The Padova- Vienna team included B. Aringer, T. Lebzelter, and 
W. Nowotny. This analysis used the COMARCS hydrostatic at- 
mosphere models, which are a further development of the origi- 
nal MARCS code of Gustafsson et al. ( 1975[ l and J0rgensen et al. 
( |1992| l. Opacities are treated using opacity tables created in ad- 
vance with the opacity generation code comaOS (Aringer et al. 
200911 ■ Compared to the MARCS models of [Gustafsson et al. 
( |2008| l, the opacities are computed with different interpolation 
schemes. The temperature and, accordingly, the pressure strat- 
ification of the models are derived under the assumption of a 
spherical configuration in hydrostatic and local thermodynamic 
equilibrium (LTE). A detailed description of the spectral synthe- 
sis is given in Nowotny et al. (2010). Sources for spectral line 



Table 4. Overview of model grid parameters used by different 
groups in Experiment 1 . 



ID" 




log 9 


[Fe/H] " 




(K)' 


(cm s^^)* 




Ml 


3600 to 40007 50 


0.5, 1.0, 1.5 


0.0, -0.3, -0.5 


M2 


3500 to 4100''/ 100 


0.5 to 2.0/0.5 


-0.5 to +0.5 / 0.25 


MS" 


2500 to 8000 


0.0 to 5.0 


-5.0 to +1.0 


M4 


3800 to 4300/ 50 


0.0 to 3.0/ 0.5 


0.0, -0.25 


M5 


2500 to 8000 / 250^ 


-1 to 5/ 0.5 


-5 to 1 / » 


M6 


3750 to 4250/ 200 


0.0 to 2.0/ 1.0 


-1.0 to 0.5/ 0.1 


Al 


3500 to 13000/250 


0.0 to 5.0/ 0.5 


-1.5 to +0.5/ 0.5 


A2 


3500 to 6000 / 250 


0.0 to 5.0/ 0.5 


-1.5 to +0.5/ 0.5 


A3 


3000 to 8000 / 100 


0.0 to 3.0/ 0.1 


0.0^ 


A4 


3500 to 6250 / 250 


0.0 to 5.0/ 0.5 


+0.5, +0.2, +0.0; 








-0.5 to -4.5 / 0.5'' 


A5 


3500 to 6000 / 250 


0.0 to 5.0/ 0.5 


0.0,-1.0 


P 


3500 to 4500/ 125 


1.0 to 2.5/0.5 


0.0, -0.5 


T 


2800 to 4000 / 100 


-0.52 to 1.34/3 


0.0 


C 


2600 to 4000 / 200 


-0.5, 0.0 


-0.5, 0.0 



Notes. For the group IDs see the subsections of Section]3] A range 
of parameters is given in the format: minimum value to maximum value 
I step size. Comma-separated values indicate a discrete set of parameter 
values. T^ff range chosen on the basis of an educated guess from the 
strength of the TiO features in the spectra, 81 models. T^ff range 
centered on T^ff from photometric calibrations. One value of [ff/Fe] 
and 2905 models for M3a, five values of [a/Fe] and 16783 models for 
M3b. Step size is lOOK between 2500 and 4000 K. In addition, a 
small set of C-rich models was used. Variable step size. [Fe/H] 
-1.5 to 0.2 / 0.1 for a subgrid with T^s 3000 to 4000 K. One model 
with a solar abundance pattern was calculated, and one model with a 
elements enhanced by +0.4. 



Recent applications to the fitting of observed high-resolution 
spectra can be found in'Lebzelter et al.'('2008') and Lederer et al. 



data are given in |Lederer & Aringer| ( |2009[ ). 



(2009 ). Here the same approach of finding the best fit was used, 
which will be briefly sketched in the following. Note that only 
the two spectra in the visual range. Star 1 and Star 2 were fitted. 
The same codes were used to calculate the two artificial sample 
spectra. Star 3 and Star 4 , and the corresponding broad band 
photometry. 

The Ml analysis started by calculating a grid of model spec- 
tra (see Table]?]) at R=300 000, rebinned to /?=80 000 afterwards. 
Broad band model colours were determined, too, following de- 
scriptions in |Nowotny e t al.'(^Ol 1 ). Microturbulence and macro- 
turbulence were set to fixed values (see Tables ]T2] and ]T3]l. The 
value for the macroturbulence was checked with the typical line 
widths in the observed spectrum, and was found to be reasonably 
chosen. 

The first step was to identify features in the optical range that 
are primarily sensitive to changes of a single stellar parameter. 
The starting point for finding the best fit was the calcium triplet 
which constrained the value for log g. This estimate was vali- 
dated by checking a few obviously log g sensitive lines between 
8500 and 8800 A. Because the strength of these lines is also de- 
pendent on the metallicity, it was decided to determine the best 
fitting log g value for all three [Fe/H] values of the grid indepen- 
dently. Next, the temperature was fixed based on the strength of 
the TiO band heads at 6652, 6681, and 6698 A. The temperature 
was derived for each of the three metallicity / log g pairs from the 
previous step. A semiempirical approach combining analysis 
with a check of the fit by eye was applied. In cases where the 
best fit seemed to fall between two points of the model grid, the 
final fit parameters were determined by interpolation between 
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these grid points. For Star 1 this resulted in the best fitting pa- 
rameter combinations (reff, \ogg, [Fe/H]) = (3900 K, 1.25, 0.0), 
(3800 K, 0.75, -0.3), and (3750 K, 0.5, -0.5). For Star 2, the 
corresponding values are (3800 K, 1.0, 0.0), (3700 K, 0.5, -0.3), 
and (3675 K, 0.5, -0.5). 

The best combination of metallicity, temperature and log g 
should then be determined by testing the overall fit of the atomic 
lines and molecular bands at shorter wavelengths. However, no 
clear decision could be made in that respect, either for Star 1 or 
for Star 2. It turned out that the fit of the spectrum below 6000 A 
became rather poor with several missing lines. Furthermore, sev- 
eral lines came out too strong from the model, and it is sus- 
pected that incorrect line blends are responsible for this prob- 
lem. Tentatively, an [Fe/H] value below solar seems to be more 
appropriate for both stars. 

To further constrain the metallicity of the star, synthetic 
colours for the chosen best fit combinations were compared with 
the observed values. For Star 1, the three given colours (see 
Table [T} are best represented by the solution for [Fe/H] = -0.3. 
For Star 2, the solution for [Fe/H] - -0.5 fits best, and the solar 
abundance model can be excluded. These solutions are given in 
Tables [T2]and[T3](see §4 below). The error bars are derived from 
the different metallicities considered. 



3.1.2. M2 

In the M2 analysis. Star 1 and Star 2 were studied by B. Plez. 
Spectra were com puted using the TURBOSPECTRUM code 
( Alvarez & Plez|19 98 ). Model atmospheres were extracted from 
the MARCS database (Gustafss on et al.| |2008). Refinement in 
the grid of models was obt ained by interpolation using the 



routines by 



Masseron 



( |2006| , provid ed on the MARCS sit^ 
Molecular line lists are described in [Gustafsson et al.| ( |2008| l, 
and atomic lines were extracted from the VALD database. As in 
the MARCS calculations, collisional line broadening is treated 
following Barklem et al. ( 2000 1, with particular broadening co- 
efficients for many Unes. 

The comparison between observations and calculations was 
done using a mean square differences computation. As a first 
step, the photometry provided with the data was used to derive a 
first estimate of T^s. Colour - T^s calibrations based on MARCS 
models were used ( |Bessell et al.||1998{ [van Eck et al.|[20TT] l. 
Spectra were then computed for spherical models of one solar 
mass, with T^g +200 K around these values, gravities typical for 
giants, a range of metallicities, and fixed values for the micro- 
turbulence parameter (see Table|4|. The differences between ob- 
served (o,) and calculated (c,) spectra with wavelength points 

were characterized by cr = computed for wide por- 

tions of the spectra (see Fig.[T|. 

This computation quickly showed that models with different 
parameter combinations (e.g. a lower T^s and a lower metallic- 
ity) could give similarly small values of cr. Also, different spec- 
tral regions would lead to best fits for different values of the stel- 
lar parameters. The reason is that the spectra are dominated by 
numerous atomic and molecular lines, for which the line position 
and strength data might not be sufficiently accurate. The result 
is then that many small differences add up to a relatively large 
cr that becomes quite insensitive to local improvements on a few 
"good" lines. Another problem is caused by the normalization of 
the observed spectra, which leaves residual slopes relative to the 
calculations. This is especially problematic for spectra of cooler 
objects with many molecular bands. Robust methods indeed use 



only the information from selected spectral ranges, where line 
data are carefully calibrated, and where spectra can be renormal- 
ized easily. This selection was not done owing to the limitations 
of present experiment. 

Values of cr were also computed for a number of macrotur- 
bulence parameters, but they tended to decrease systematically 
with increasing macroturbulence. This is again explained by the 
presence of numerous lines, many of which were not well fit- 
ted. It proved more efficient to use values derived from an eye- 
inspection of the spectra (see Tables 12 and 13 i. The best set of 
parameters (Teft, log g, metallicity) was derived for the various 
spectral regions (Fig.[T]i from the model giving the smallest cr, 
and for the few models giving similarly small values. Later, an 
inspection by eye assigned higher weight to regions that gave a 
better global fit (and a correspondingly smaller cr). The results 
for Star 1 and Star 2 are given in Tables 12 and 13 and are dis- 
cussed below. 



Star 1 and Star 2: Photometry suggested the T^s values given 
in Tables [12] and [13] The values were consistent for J - K and 
V -K in the case of Star 2, but not for Star 1. The preferred cali- 
bration for cool stars is Ti^fi-{V -K), which is quite insensitive to 
metallicity, and is strongly dependent on T^g. It is, however, also 
sensitive to the C/O ratio, especially when C/O>0.9; this is the 
reason for using a combination ofV — K and J — KXo determine 
S-type star parameters ( van Eck et al.|20lT] l. 

For Star 1, the best cr- values for the first four spectral win- 
dows (Fig.[l]| were 0.12, 0.047, 0.055, and 0.055. The number of 
models with a cr within 10% of the smallest were 35, 5, 13, and 
7. This procedure yielded an estimate of Tss of about 3800 K and 
log g around 1.5, but there were other models with Tetf =3700 K 
to 3900 K, and log g from 0.5 to 2.0, which also gave small val- 
ues of cr . The metallicity was not well constrained (solar, or 
maybe sub-solar). Attempts were made for the Ca II lines (8500- 
8700 A) in the Gaia Radial Velocity Spectrometer (RVS) range, 
leading to a best fit for reff=3850 K, log g-1.25, and solar metal- 
licity. The individual Ca abundance was also varied, once the 
other model parameters were fixed, and was found to be solar 
to within 0. 1 dex. The fit of the 5000-5200 A region (Mg I and 
MgH lines) gave consistent stellar parameters, and a solar Mg 
abundance. For Star 2, a similar approach resulted in the param- 
eters given in Table 
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This gave a best fit in the 6400 A region 
but not in the Ca II IR triplet region, unless Ca is shghtly over- 
abundant (h-0.2 dex). 

Star 3 and Star 4: Photometry gave the Tefr values given in 
Tables [T4l and [TS] No further work was done on these stars. The 
Teff values were consistent for J - K and V - K in the case of 
Star 3, but not for Star 4. A consistent Teff =3400 K could have 
been obtained from both colours, if the star is s-element rich 
(about -1-2 dex relative to the Sun). As this toy-star had a solar 
composition, this points to a difference between MARCS and 
COMARCS model spectra, presumably due to differences in the 
opacities. 

3.1.3. M3 

The M3 team consisted of C. Worley, P. de Laverny, A. Recio- 
Blanco and G. Kordopatis. In this case, two separate spectral 
analyses were carried out for each of Star 1 and Star 2 us- 
ing pre-existing procedures that have been established for two 
separate research projects. Both projects consist of automated 
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pipelines that feed the spectra into the stellar classification al- 
gorithm MATISSE (MaTrix Inversion for Spectrum SynthEsis). 
Thus, the two stars were analyzed in a completely blind way us- 
ing automated procedures. The spectra were degraded to much 
lower resolutions (R=s6500 for M3a and R^15000 for M3b) 
and no photometry was used to obtain prior estimates of the 
stellar parameters. Only the observed spectra were used in the 
search across a wide range of stellar parameters (dwarf to gi- 
ant, metal-rich to metal-poor FGKM stars). The entire process, 
during which the spectra undergo wavelength selection, cosmic 
ray cleaning, radial velocity determination and correction, nor- 
malisation, and analysis in MATISSE, lasts only a few minutes 
for each pipeline. A key feature is the iteration between the nor- 
malisation procedure and the stellar parameter determination in 
MATISSE, whereby synthetic spectra generated for the preced- 
ing set of stellar parameters are used to normalise the observed 
spectrum for the next parameter determination. In this manner 
there is convergence to the final set of stellar parameters, the cor- 
responding synthetic spectra, and the final normalised observed 
spectra. 

The MATISSE algorithm, initially developed to be used in 
the analysis of the Gaia RVS spectra, is based on a local multi- 
linear regression method ( Recio-Blanco et al. 2006 Bijaoui et al 
2008|. A stellar parameter 6 (in this case: T^ff, log g, or [M/H]) of 
a star is determined by the projection of the observed spectrum, 
0(A), onto a vector function Bg{A). The Bg(A) vector is an opti- 
mal linear combination of spectra, S i{A), in a grid of theoretical 
spectra. The product - jj 2 Bg{A) 0(A) is calculated. 

A=1,N 

In the training phase of MATISSE the Bg(A) vectors are cre- 
ated from the grid of synthetic spectra, with Bfl(/l) - 2 Oj Si(A), 
and a, being the weight associated with the spectrum S , giving 
the maximum correlation between 0, and 0, in the training grid. 
The sensitivity of a wavelength region to the particular stellar pa- 
rameter 6 is reflected in the corresponding Bg(A) vector. The syn- 
thetic spectra cover the entire optical domain and have been built 
using flie spectral line format ion code TURB OSPECTRUM 
(Plez, private communication; Alvarez & Plez 1998|) a nd the 
MARCS stellai- atmosphere models (|Gustafsson et al.||2008)l. 
We decided not to perform any convolution (spectral resolution 
and instrumental profile, stellar rotation, macro turbulence,...) for 
these spectra in order to keep an easy and fast adaptation to the 
properties of any spectrograph and/or astrophysical application 
(see also |de Laverny et al. 2012) . At the comparably low spectral 
resolution the M3 analyses were performed, these broadening 
parameters were of negligible relevance. 

This grid covers the range of parameters given in Table |4] 
The atomic line lists were taken from VALD (August 2009) and 
for M3a were calibrated to the Sun and Arcturus. The molec- 
ular line lists were provided by B. Plez and included CH, OH, 
MgH plus several isotopic compositions, SiH, CaH, FeH, C2, 
CN, TiO, VO, and ZrO. The grid is presented in more detail in 
de Laverny et al. (2012, to be submitted). For each of the two 
research projects there was a training phase during which the 
Bg(A) functions for each spectral set-up (wavelength regions and 
resolution) were generated from the initial synthetic spectra grid. 



FLAMES/GIRAFFE MATISSE Analysis: The first analysis 
method used (M3a) has been developed for the analysis of 
~700 galactic disk stars o bserved using FLAMES/GIRAFFE 
(Kordopatis et al.||2011a|b)l. The wavelength domain and reso- 
lution of this project correspond to that of the Gaia RVS (8470- 
8740 A and R~6500), of which a key spectral feature is the Call 



Infra-Red Triplet (LR8 setup of GIRAFFE). For this analysis 
some small features were masked in order to eliminate sky lines 
and the cores of some strong lines which were not well repro- 
duced by the synthetic spectra. Also the line list was calibrated 
to the Sun and Arcturus. For this analysis standard galactic a el- 
ement enhancements were assumed for the stellar atmospheric 
models; this is the reason for the lower number of nodes in the 
synthetic spectra grid, compared to M3b (see below). For a full 
description of the FLAMES/GIRAFFE analysis including the 
adaption of the line lists and the error analysis, see Kordopatis 
|etar] ( |2011a|bl l. 

AMBRE-FEROS MATISSE Analysis: The second approach 
(M3b) has been designed for the analysis of the archiv ed spec- 
tra of the FEROS spectrograph (AMBRE ProjecQ; see Worley 
|etaL] ( |20T2| for details. 

The grid of synthetic spectra for the AMBRE Project in- 
cludes the full range of parameters, and the spectra cover the 
entire optical domain (de L averny et al.|2012| l. Due to the large 
wavelength range no calibration of the line list was carried out. 
Comprehensive comparisons were made between the AMBRE- 
FEROS stellar parameters and corresponding stars in the S^N li- 
brary, the PArametres STELlaires (PASTEL) database ( Soubiran 
et al.|2010 ) and the dwarf stellar sample in Bensby et al. ( 2003| l. 
For the combined sample of 178 stars from S'^N and PASTEL, 
the dispersions in Teff, log g, and [M/H], are 150 K, 0.25 dex. 



and 0.12 dex respectively. For the 66 stars in Bensby et al. 
( 2003| l the dispersions in Teff, logg', [M/H], and [a/Fe] are 86 K, 
0.17 dex, 0.06 dex and 0.04 dex respectively. Hence there is 
excellent agreement between AMBRE-FEROS and these high 
quality, non-automated, spectroscopic analyses. 

For the AMBRE-FEROS analysis of Star 1 and Star 2 
(M3b), as the spectra did not encompass the full AMBRE- 
FEROS wavelength range, only the regions between 4900- 
6730A were extracted (~435A in total). The spectra were then 
convolved to R~15 000. High resolution Hinkle atlases of the 
Sun and Arcturus ( [Wallace et al.| 2007; Hi nkle et al.[ 2000) were 
analyzed concurrently with Star 1 and Star 2. The parameters 
that were determined are T^s = 5783±80 K, loggr = 4.30+0.15, 
[M/H] = O.OliO.lO, [a/Fe] = 0.00+0.10 for the Sun, and 
r^ff = 4306+80 K, log^ = 1.80+0.15, [M/H] = -0.66+0.10, 
[of/Fe] - 0.24+0.10 for Arcturus, in excellent agreement with 
accepted values. 

Results for Star 1 and Star 2 for both approaches are listed in 
Tables [T2]and[T3] Both analyses confirm that these stars are cool 
giants. For Star 1 and Star 2 both the M3a and M3b analyses 
show very good agreement in T^s and [Fe/H] in comparison to 
the stellar parameters quoted in Table [T] For Star 1 there is an 
absolute difference in gravity of Alogg = 0.60 between M3a 
and M3b. However both values differ by only 0.3 dex from the 
adopted value stated in Table [T] 

For the M3a analysis the small spectral domain of 8470- 
8740A at this fairly low spectral resolution is known to con- 
tain rather poor (and few) spectral signatures that are sensitive 
to gravity, leading to degeneracies in the derivation of the stel- 
lar parameters (see [Kordopatis et al.|[201 la| l. Gravity problems 
for the Ca II IR triplet arise specifically for cool dwarfs (lack of 
signatures) and subgiants (Teff-log g deceneracy). The problem 
is much less important for the giant stars studied in this paper. 



' performed under a contract between the Observatoire de la Cote 
d'Azur (OCA) and the European Southern Observatory (ESO). 
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For the M3b analysis it was determined that the normalisa- 
tion of the magnesium triplet at ~5170A was driving the issues 
with the gravity determination. At these cool temperatures this is 
a difficult region to normalise and, as the Mgib triplet is a domi- 
nant contributor to the gravity determination, the less than ideal 
normalisation of this region at each normalisation/MATISSE it- 
eration drove the automated process to less precise gravity val- 
ues, in particular for Star 2. In the full AMBRE-FEROS analysis 
a three times wider wavelength interval is used which includes 
significant, and better distributed, contributors to the gravity. 
This reduces the impact of the MgIb triplet at these cool tem- 
peratures. 

3.1.4. M4 

The M4 analysis was done for Star 1 by K. Eriksson. A grid 
of model atmospheres was calculated for the parameters given 
in Table |4] All models were spherical with one solar mass. For 
more details on the MARCS model atmospheres seejGustafsson 
et al.| ([2008 ). For each of these model atmospheres synthetic 
spectra were calculated, for the three wavelength regions indi- 
cated in Fig. [l] with a wavelength step of 0.02 A using the line 
formation code BSYN. The spectra were then downgraded to a 
resolution of 80000. 

For the atomic lines the VALD data base was used. For com- 
puting the atmosphere models, absorption data for transitions of 
18 molecules were included (see Table 2 of [Gustafsson et aT] 
|2008 ). The molecular species included in the synthetic spectrum 
calculations are listed in Table 2] 

Table 5. Molecular species used in the M4 synthetic spectrum 
calculations. 



Species 4900-5400 A 


6100-6800 A 


8400-8900 A 


C, • 






CN . 


• 


• 


MgH . 
SiH • 


• 




TiO • 


• 


• 


CaH 


• 




ZrO 


• 




FeH 




• 



with different curves for diff'erent logg values. The best fit was 
obtained for an effective temperature of 4000 K or slightly less 
and for a loggr of 3.0 or less. The same procedure for model 
atmospheres with a metal abundance of [Fe/H]= -0.25 yielded 
the result shown in Fig. |2] (bottom left). Note that the best fit is 
shifted to lower T^s values and lower log g values by 100 K and 
1 dex, respectively. 

8400-8900 A: In this wavelength region, very interesting from 
a Gaia point of view, we find the Call infrared triplet, which 
is strong in cool giant stars. The plots for the 8475-8575 A 
wavelength interval for solar abundances and for [Fe/H]= -0.25 
can be found in Fig.|2](top and bottom center, respectively). The 
best fit for solar abundances is shown in Fig. [3] This region is not 
very sensitive to the effective temperature, and quite sensitive to 
the surface gravity. Using a lower metal abundance resulted in a 
lower surface gravity. 



+0.2 




I ' ' ' ' ' 1 

8480 8490 8500 8510 8520 

Wavelength [A] 

Fig. 3. M4 analysis of Star 1: The best fit for the 8475-8575 A 
region, showing a Call IR triplet line, for solar abundances 
(refF=4050K, log gf=1.5). Bottom: Scaled flux, observed and cal- 
culated; top: observed minus calculated flux. 



The flux was adjusted by a continuum scale factor deter- 
mined from 100 A- wide regions for the different parts of the 
spectrum. The synthetic spectra were then cross-correlated with 
the given observed spectrum for Star 1 for a number of wave- 
length regions and a;t'^-value was calculated {x^ = YjhA , 
where AA is the wavelength interval, m is the model spectrum, 
and o is the observed spectrum). The minimum ;t'^-value should 
then give the best model atmosphere for that wavelength region 
if the included line lists are complete and accurate enough. The 
results for some representative wavelength regions are as fol- 
lows. 



6100-6800 A: The region 6400-6700 A was used to set the 
best flux scaling of the observed spectrum. The scale factor de- 
pended to some degree on the model atmosphere used, but did 
not vary by more than a few percentage units. With this scal- 
ing factor the wavelength region 6100-6400 A was investigated. 
The result for solar abundances can be seen in the top left panel 
of Fig. [2] where the ;i'^-values are presented as a function of Teff 



4900-5400 A: This wavelength region has a much weaker 
temperature- and gravity sensitivity, as can be seen in Fig. |2] 
(right) for the interval 5150-5200 A containing the MgIb triplet 
lines. The "best" fit for the [Fe/H]= -0.25 models is shown in 

Fig.g 

In conclusion, within the M4 analysis, the 6100-6400 A re- 
gion, with its temperature-sensitive TiO bands, was best suited 
for estimating the effective temperature for Star 1. The 8500 A 
region was best suited for estimating the surface gravity. The 
solutions for the two adopted metallicity values are given in 
Table 12 It should be noted that the M4 analysis did not result 
in a unique set of best-fit parameters: the temperature-sensitive 
region gave a higher log g than the gravity-sensitive region, and 
vice versa. The blue-green region (5150-5200 A) gave both a 
higher Teff and a higher log g than the other two regions. Most 
of the other tested wavelength intervals behaved similarly to the 
blue-green region, i.e. showing very broad minima at rather high 
Teif and loggr in the plots. These were four 50-100 A wide 
intervals between 4900 and 5400 A, and a similar number of in- 
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4100 4200 4300 



3800 3900 4000 4100 4200 4300 




3800 3900 



4100 4200 4300 



Fig. 2. M4 analysis of Star 1: values as a function of T^s for different log g values (see legend) and for three different wavelength 
regions (left, middle, right). Top row: solar abundances, bottom row: [Fe/H]= -0.25. 




5150 



5160 



5170 5180 
Wavelength [A] 



5190 



5200 



Fig. 4. M4 analysis of Star 1: The best fit for the 5150-5200 A 
region, showing the Mg lb triplet lines, for sub-solar abundances 
([Fe/H]= -0.25; reff=4150 K, logg=3.0). Bottom: Scaled flux, 
observed and calculated; top: observed minus calculated flux. 



tervals in the region 8400-8900 A, which did not contain any of 
the Call triplet lines. 



An analysis of individual (e.g. Fel) lines to derive an esti- 
mate of [Fe/H] was not attempted. Neither was a analysis in 
the metallicity dimension (although one can note that these val- 
ues are on average somewhat smaller for the -0.25 models in 
Fig.|2]i. A further conclusion is that this procedure is to some de- 
gree dependent on the scale factors, i.e. it is important to have a 
good fit to the stellar continuum (which is hard to do, due to the 
effect of molecular bands). 

3.1.5. M5 

This part of the experiment was led by C. Abia. The M5 anal- 
ysis focused on the coolest stars of the experiment. Star 3 and 
Star 4, because the M5 team members are more familiar with 
the spectral region to be analyzed in these stars (the //-band) 
where, furthermore, telluric lines are absent. 

The original atomic line list was taken from the VALD 
database (version 2009). For the molecules, C2 lines are from 
WahHn & Plez| (pOOS]), CO fines come from |Goorvitchl ( [1994) 1, 
an the CN and CH lines were assembled from the best avail- 
able data as described in Hill et al. ( 2002 1 and |Cayrel et al 
(2004i. The molecular lists also include lines of OH, TiO, CaH, 



SiH, FeH and H2O taken from the HITRAN database (Rothman) 
|etal.|2005) . First, the atomic and molecular lines were calibrated 
by obtaining astrophysical g f values using the high resolution 
Solar ( [Livingston & Wallace|199T| l and Arcturus ( [Hinkle et al.| 
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|1995| l spectra in the 1.5/im region. In both analyses MARCS 
atmosphere models ( [Gustafsson et al.|20081 ) were used, with pa- 
rameters reff/logg/[Fe/H]=5777/4.4/0.0 and 4300/1 .5/-0.5 for 
the Sun and Arcturus, respectively. The reference abundances 
were th ose by|Asplund et al.| ( 2009) l for the Sun, and those de- 
rived by Peterson et al.| ( 1993| for Arcturus. The 1 .5yum region is 
dominated by CO, OH, CN and C2 absorptions (in order of de- 
creasing importance). The main features were calibrated on the 
Arcturus spectrum. Synthetic spectra in LTE were computed us- 
ing the TURB OSPECTRUM V9.02 code described in |Alvarez 
Plez ( 1998[ l. The theoretical spectrum was convolved with a 
Gaussian function with a FWHM of ~300 mA. The stellar pa- 
rameters were estimated from the available photometry using the 
recent photometric calibrations by Worthey & Lee ( 2011[ l (see 
TablesfT4landfT5ll. 



Star 3: A spherical (1 M©) model atmosphere with parameters 
determined from the photometric colours was taken from the 
grid of [Gustafsson et aL] (p008). The fit (by eye) with this set 
of parameters to the observed spectrum was quite good in the 
full spectral range. The simultaneous fit to some OH, CO and 
CN lines allowed an estimate of the CNO abundances, namely 
C/N/0= 8.37/7.90/8.70, i.e. C/0= 0.46. A rough estimate of 
the metallicity from fits to some metallic lines in the region was 
compatible with [Fe/H] ^ 0.0 + 0.2. Finally, according to the 
line list used, some CO lines are sensitive to variations in the 
'^C/'^C ratio. Several fits unambiguously resulted in a carbon 
isotopic ratio lower than the solar value (89). However, due to 
the weakness of the '^CO lines only a lower limit was derived, 
i2c/i3C> 25. 



Star 4: With a model atmosphere with parameters determined 
from the photometric colours it was impossible to fit the ob- 
served spectrum. In particular, the predicted intensity of the CO 
lines was lower than observed. In this case, the stellar parameters 
were estimated through a test comparing observed and theo- 
retical spectra computed for several choices of Teff/log g/[Fe/H] 

t^^(') = ^ c@+o\'o ' where c is the computed spectrum, and o is 
the observed spectrum). 

The best fit was obtained with a spherical 2 M© MARCS 
model (K. Eriksson, private communication) of parameters 
7'eif/log0/[Fe/H]= 3500/0.0/0.0 (x^ < 1.5 over 20-30 A re- 
gions). With these stellar parameters the corresponding photo- 
metric colours according to [Worthey & Lee] ( |2011 1 are instead 
(y - /) = 2.49, (J-K)^ 1.18, and (V - K) ^ 5.15. In the next 
step, the CNO abundances were estimated in a similar way to 
Star 3, resulting in C/N/O = 8.52/7.88/8.68, i.e. C/0= 0.69, 
and a lower limit of 35 for the '^C/'^C. The estimated metallic- 
ity was also compatible with [Fe/H]= 0.0. The microturbulence 
parameter was set to a fixed value (see Table [T5|. Note that the 
C/O ratio in the model atmosphere used (close to 1) was not 
consistent with the C/O ratio derived. The fact that this M-type 
star has some carbon enrichment may explain why a theoretical 
spectrum computed with the stellar parameters giving the best fit 
to the photometric colours does not fit the observed spectrum. In 
cool carbon enhanced stars, the actual C/O ratio is a critical pa- 
rameter which determines the structure of the atmosphere affect- 
ing also the photometric colours. The results of the M5 analysis 
for both stars are summarized in Table [TSl 



Table 6. Iron line selection for the M6 stellar parameter deter- 
mination (Sect. 3.1.6| l. For each line, the wavelength A, the lower 
level excitation potential Eio^, the gf value, and the source for 
the gf value are given. 
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K07 
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00 i J i 




— 1 70^ 


JVU / 


881.689 


4.988 


-2.203 


K07 


883.402 


4.218 


-2.558 


K07 


884.674 


5.010 


-0.111 


K07 


889.140 


5.334 


-1.108 


K07 


Fell 








611.332 


3.221 


-4.230 


RU 


614.774 


3.889 


-2.827 


RU 


641.692 


3.892 


-2.877 


RU 


643.268 


10.930 


-1.236 


RU 


645.638 


3.903 


-2.185 


RU 


651.608 


2.891 


-3.432 


RU 



Notes. K07 ...|Kurucz[ 
[Raassen & Uylings| ( |1998y 



j2007f , BWL ... [O'Brian et"ar] ( |l99l| ), RU . 



3.1.6. M6 

The M6 team, consisting of T. Merle and F. Thevenin, did not 
try to reproduce molecular lines and chosen to work with only 
Star 1 which has less affected spectra. For the M6 analysis, the 
photometry information was used to get starting values for Tetf, 
loggf, and [Fe/H]. With the proposed approximate colours and 
using the theoretical calibrations by [Houdashelt et aL ( 2000[ l, 
the Teff and log g values given in Table|12|were obtained. For K 
giant stars, a good approximation is to take a mass of 1 Mq. The 
spectroscopic analysis started with the selection of Fe I and Fe II 
lines. Care was taken to remove strong lines, multiplets with no 
dominant log gf (> 10%) and polluted or unclear line shapes 
in the observed spectrum. 21 Fel lines were chosen in the near- 
IR range. Unfortunately, all Fe II lines in the near-IR part were 
too weak to use for a fit. Then, six Fell lines were selected 
in the [6100 - 6800]A range. For each line, radiative damping 
was estimated using the Unsold formula; hydrogen elastic col- 
lisional damping came from ABO theory (Anstee & O'Mara 
1995; Barkle m & 0'Mara|[T997l [Barklem et a l. 1998) when 
available (e.g. the six Fell lines), otherwise, classical Unsold 
theory was used with an enhancement factor of 1.5 (e.g. the 21 
Fe I lines). The line list is given in Table|6] 

The final values of the stellar parameters for the M6 anal- 
ysis were constrained by the method of ionization equilibria 
between Fel and Fell and excitation equilibria for Fel. For 
this analysis, the MOOG2009 code ( |Sneden||197"3] l with spher- 
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ical MARCS model atmospheres ( |Gustafsson et al.||2008| l and 
atomic line lists from VALD were used. The inconsistency be- 
tween spherical atmosphere and plane parallel radiative transfer 
is negligible (Heiter & Eriksson |2006| l. The chemical compo- 
sition was from iGrevesse et al.| ( 2007| . For each model atmo- 
sphere the abundance was determined for each line by fitting the 
observed spectrum while varying the iron abundance in steps of 
0.1 dex around the abundance value adopted for the model at- 
mosphere. The continuum level was determined using a local 
scaling factor Then, [Fel/H] was plotted as a function of the 
lower level excitation potential £, of the lines, and using linear 
regression, the model with the smallest slope was selected. For 
the ionization equilibrium, the mean [Fe I/H] and [Fe II/H] was 
compared for each model, and the model with the smallest abun- 
dance difference was selected. The best model atmosphere had 
parameters as given in Table 12 The uncertainties come from 
the large steps in the stellar parameters (200 K for red , 1 for 
log g , and 0.1 dex in metallicity) taken to perform the analysis. 
With these large steps, the same atmospheric parameters were 
obtained both from the slope of [Fe/H] vs £, and from the ion- 
ization equilibrium. 

The radial velocity was determined using 21 weak and rela- 
tively unblended lines in the [6100 - 6800] and [8400 - 8900] A 
domains. The lines were synthesized using the MOOG2009 LTE 
code and shifted in steps of 0.5 km s ' for all the lines selected. 

3.2. ATLAS model atmospheres (A) 
3.2.1. Al 

This working group consisted of G. Wahlgren and R. Norris. 
Using synthetic photometric colours from Kucinskas et al. 
( 2005| l, initial estimates of temperature, gravity, and metallicity 
were made. With a least sum squares approach, the colours of 
the provided stars were compared to the colours in Kucinskas' 
table, and the parameters of the best matching synthetic colours 
were interpolated to the given stars (see Tables [T2| to [TSjl. 

Using these estimates of stellar parameters, appropriate 
ATLAS9 model atmospheres were selected from the grid of 



Castelli & Kurucz (2003 , and synthetic spectra were created 



with SYNTHE ( |Kurucz| |1993b| . ATLAS is an LTE, plane- 
parallel model atmosphere code. Molecules from the Kurucz 
suite of line lists d Kurucz & B"ein995 hydrides, CN, CO, TiO, 
SiO, H2O) as well as other molecules were included in the sta- 
tistical equilibrium calculations. The use of the gravities and 
metallicities obtained from broadband colours produced syn- 
thetic spectra which did not match the workshop spectra (at lines 
known to be sensitive to gravity and metallicity). Therefore, spe- 
cific spectral lines as indicators of gravity and metallicity were 
used. The Al analysis was continued for the K-type stars Star 1 
and Star 3, because the contributing authors were not convinced 
that they had sufficient molecular opacities for late M-type stars. 

Atomic line data in this analysis come from Kurucz (approx- 
imately 2008). Tables |7]and[8]list the lines used in the abundance 
determination for Star 1 and Star 3, respectively. Lines from neu- 
tral species only were used. 



Star 1: Synthetic colours suggested a cool, metal deficient gi- 
ant or supergiant. Comparison to UVES spectra showed close 
matches with both HR 7971 (K3 II/III) and HR 611, which 
has ben classified as M0.8 III from nan^ow-band TiO/CN pho- 
tometiy (Wing 1978). The infrared Can triplet (8498, 8542, 
8662 A) and the Fe i calibration Unes at 8327, 8468, 8514, and 



Table 7. Lines used in the Al analysis to determine abundances 
for Star 1. For each line, the atomic number Z, the wavelength 
A, the gf value, the lower level excitation potential E\o^, and the 
line abundance are given. 



z 


A 


log(R^) 


Source" 




Abundance 




(nm) 






[cm-'] 


(log(N„)= 12.00) 


14 


568.448 


-1.650 


GARZ 


39955.053 


6.95 


14 


569.043 


-1.870 


GARZ 


39760.285 


7.10 


14 


570.110 


-2.050 


GARZ 


39760.285 


7.08 


14 


577.215 


-1.750 


GARZ 


40991.884 


7.04 


14 


579.307 


-2.060 


GARZ 


39760.285 


7.15 


20 


558.875 


+0.210 


NBS 


20371.000 


4.83 


20 


559.011 


-0.710 


NBS 


20335.360 


5.23 


20 


559.446 


-0.050 


NBS 


20349.260 


4.81 


20 


560.128 


-0.690 


NBS 


20371.000 


5.13 


20 


615.602 


-2.200 


NBS 


20335.360 


5.88 


20 


645.560 


-1.350 


NBS 


20349.260 


5.48 


20 


817.329 


-0.546 


K88 


39349.080 


5.90 


20 


820.375 


-0.848 


K88 


36575.119 


5.90 


20 


829.218 


-0.391 


K88 


39340.080 


5.80 


22 


546.047 


-2.880 


SK 


386.874 


4.55 


22 


556.271 


-2.870 


MFW 


7255.369 


4.55 


22 


802.484 


-1.140 


MFW 


15156.787 


4.45 


22 


806.825 


-1.280 


MFW 


15108.121 


4.65 


22 


841.229 


-1.483 


MFW 


18482.860 


4.55 


22 


867.527 


-1.669 


MFW 


8602.340 


4.50 


22 


868.299 


-1.941 


MFW 


8492.421 


4.50 


22 


869.233 


-2.295 


MFW 


8436.618 


4.70 


22 


873.471 


-2.384 


MFW 


8492.421 


4.70 


26 


552.554 


-1.330 


FMW 


34121.580 


6.70 


26 


552.890 


-2.020 


FMW 


36079.366 


7.30 


26 


552.916 


-2.730 


FMW 


29371.811 


7.10 


26 


553.275 


-2.150 


FMW 


28819.946 


6.70 


26 


554.394 


-1.140 


FMW 


34017.098 


6.70 


26 


555.798 


-1.280 


FMW 


36079.366 


6.65 


26 


556.021 


-1.190 


FMW 


35767.561 


6.75 


26 


829.351 


-2.126 


K88 


38678.032 


6.90 


26 


832.705 


-1.525 


FMW 


17726.981 


6.70 



References as listed in [Kurucz & Bell| ( |l995| l: 



Notes. 



http://kurucz.harvard.edu/LINELISTS/LINES/gfall.refl 



8689 A noted in 



Keenan & Hynek| ( |1945| were used to deter- 

4000 K. 



mine loggf of 1.5 in a temperature range of 3900 
Deficiencies in iron (see Table [T2] ) and the a elements were 
of particular note ([Si/H]=-0.45 + 0.03, [Ca/H]=-0.90 ± 0.15, 
[Ti/H]=-0.38 ± 0.03). The lines used in these measurements are 
listed in Table |7] 

Of the a elements, the abundance of calcium resulting 
from the Al analysis is particularly low. For the six lines be- 
low 8000 A, which have a lower excitation energy close to 
20 000 cm"', weaker gf/-values correlated to higher abundances. 
For example, the 5589 A and 5601 A lines have the same 
lower excitation energy but have different ^/-values. Of these, 
5601 A with a log(gf/)=-0.690 dex suggests an abundance of 
5.13 whereas 5589 Awith log(gr/)=0.210 dex suggests an abun- 
dance of 4.83. This trend does not exist for the three calcium 
lines above 8000 A which were included in the abundance de- 
termination. Although each of these lines has a more energetic 
lower excitation energy, each suggests a higher abundance than 
the shorter wavelength lines with similar ^^/-values. These three 
p988) calculated gr/-values as opposed to the 
.[1969 1 values used for the shorter wavelength 



lines use Kurucz 



NBS (Wiese eta 



lines. The |Kurucz| ( |2007| l calculated values for the four lines be- 
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Table 8. Fel lines used in the Al analysis to determine iron 
abundances for Star 3. For each line, the wavelength A, the gf 
value, the lower level excitation potential fiiow, and the line abun- 
dance are given. 



A 


log(g/) 


Source" 




Abundance 


(nm) 






[cm-'] 


[log(N„)= 12.00] 


1549.034 


-4.574 


O 


17726.987 


7.85 


1558.826 


+0.323 


K94 


51359.489 


7.70 


1561.115 


-3.822 


K94 


27523.001 


7.65 


1564.851 


-0.714 


K94 


43763.977 


7.90 


1565.283 


-0.476 


K94 


50377.905 


7.80 



References as listed in 'Kurucz & Belli jl995f 



Notes. 

http://kurucz.harvard.edu/LlNELlSTS/LlNES/gfall.ref" 

low 6000 A have larger gr/-values than NBS reports. Presuming 
that this correlation between calculated and experimental values 
would continue for longer wavelength lines, this indicates that 
for the three lines over 8000 A, the reported gf/-value is an up- 
per bound. If the correlation between higher abundance and low 
gr/-value in the shorter wavelength values is not entirely the re- 
sult of poor atomic data; it is likely the result of non-LTE effects, 
to which calcium lines are particularly sensitive (for a non-LTE 
analysis of a different set of calcium lines see Section |5.3| l. The 
results of the Al analysis for Star 1 are given in Table |12| 

Star 3: Synthetic colours suggested a metal deficient giant of 
effective temperature 4300^400 K. Comparison of the work- 
shop spectrum with synthetic spectra produced with parameters 
obtained from broadband colours showed that while the \ogg 
and temperature indicated by synthetic colours fit, the metal- 
licity was closer to the solar value than the synthetic colours 
suggested. The Al team determined abundances of carbon, ni- 
trogen, and oxygen for Star 3 by fitting molecular lines. There 
are strong OH and CN lines present, as well as a CO bandhead, 
in the wavelength region of the experiment. Unblended OH lines 
were used as indicators of the oxygen abundance and CN and CO 
lines as indicators of carbon and nitrogen abundances. Several 
iron lines, listed in Table |8] served for the determination of the 
iron abundance. Despite hints from synthetic colours that the star 
was metal deficient, carbon, nitrogen, oxygen, and iron were all 
enhanced. Lines of some other elements, though not analyzed, 
suggested similarly enhanced abundances. The star was found 
to be oxygen rich. The results of the Al analysis for Star 3 are 
summarized in Table [T4l 



3.2.2. A2 

The A2 team (consisting of J. Maldonado, A. Mora, and B. 
IVIontesinos) decided to work only with the optical region of the 
spectrum, and therefore, only Star 1 and Star 2 were analyzed. 

To determine the radial velocity, the spectra of the target stars 
were cross-correlated against spectra of several radial velocity 
standards (Maldonado_et al. 20 10). Cross-correlation was per- 
formed with the IRApj^] task ^cor. Spectral ranges with promi- 
nent telluric lines were excluded from the cross-correlation. 

A grid of synthetic spectra was calculated for the param- 
eter range given in Table [4] To compute the spectra, team 

IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the National 
Science Foundation. 



A2 used ATLAS 9 model atmospheres and the SYNTHE code 
(Kurucz' 1993a"b'), adapte d to wor k under the Linux platform by 
jSbordone et al. ( 2004 ) and Sbordone (2005 Line data are from 
the Kurucz web site (2005 version). The following molecules 
are included: C2, CH, CN, CO, H2, MgH, NH, OH, SiH, SiO. 
Atomic line data are lis ted in Table [9| The new opacitiy dis- 
tribution functions from [Castelli & Kurucz] ( |2003| ) were used. 
A mixing length parameter of 1 .25 was used. All spectra were 
computed with a resolution of 300 000, 

To obtain the "best parameters" of the target stars, A2 com- 
pared the equivalent widths (EWs) of a sample of spectral lines 
measured in the target stars with the EWs of the same lines mea- 
sured in each synthetic spectrum. To compile a list of "well 
behaved" lines, the A2 team calculated synthetic spectra for 
two stars with accurately known stellar parameters (namely, the 
Sun, G2V, and Procyon, F51V) and compared them with high- 
resolution observed spectra. Relatively isolated lines, with a rea- 
sonable stretch of flat continuum around the limits, were selected 
to estimate the EWs with high accuracy. Since in all cases the 
synthetic lines fit the observed spectrum, it can be expected that 
the atomic parameters are fairly reliable. The list of lines used is 
given in Table[9] The abundance ratios of the individual elements 
relative to each other were kept fixed to the solar ratios. Final 
parameters were obtained by using a reduced fitting method. 
Lines with large EWs (>100 mA) were not used. The results are 
listed in Tables [T2l and [T3] 

The chosen method clearly can be further improved to ob- 
tain more accurate parameters. The first problem identified dur- 
ing the A2 analysis was the line selection. Although an attempt 
was made to compute a list of well behaved lines, it is clear that 
such a list depends on the spectral type of the target stars. The 
selection was started using as reference the Sun and Procyon, 
but some well behaved lines in the Sun were not well behaved 
in Procyon and vice versa. In addition, since the program stars 
of this experiment were cooler than the chosen reference stars, 
the selection had to be reviewed several times in order to avoid 
blended lines or lines not present in the program stars. 

This is related to the second problem, namely how to mea- 
sure EWs. Although EWs can be measured "by hand" (for exam- 
ple using the IRAF task splot), this is not feasible even for this 
experiment with only two target stars, since more than six hun- 
dred synthetic spectra need to be analyzed. The A2 team devel- 
oped a code which performs an integration around the center of 
each line using a fixed width. This could definitely be improved. 
Ideally, the code should be "intelligent enough" to decide on the 
width of the integration band according to the line profile. An 
alternative option, on which the A2 team is still working, is to 
compute EWs using model atmospheres and an abundance com- 
putation program such as WIDTH9 ( )Castelli|2005| ). 

3.2.3. A3 

For A3, H. Neilson attempted to determine the effective temper- 
ature, gravity, iron abundance, and carbon-to-oxygen ratio for 
Star 3 and Star 4 based on the broad-band colours and spec- 
tra provided. Model stellar atmospheres and synthetic spectra 
were computed using Fortran 90/95 versions of the ATLAS code 
( [Lester & Neil son 2008 ). The new versions of the code can com- 
pute model stellar atmospheres assuming either plane-parallel 
or spherically symmetric geometry and either opacity distribu- 
tion functions or opacity sampling. Lester & Neilson (2008| 
demonstrated that model stellar atmospheres computed with 



http://wwwuser.oat.ts.astro.it/atmos/Download.html 
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Table 9. Lines used by the A2 participants (Sect. 3.2.2 1. For each line, the species code, the wavelength A and the gf value are 
given (taken from the line lists of Kurucz). Lines included in the final fit have been marked with the symbol *. The species code 



follows Kurucz's notation, e.g. for iron lines (atomic number 26); 26.00 - Fel, 26.01 =Fe II, etc. References as listed in Kurucz & 
|Bell| ( |1995] l° . Furthermore, the table gives the measured equivalent widths of the lines in the two benchmark stars. 



Ion 






kj \J Li J. ^ 


FW Star 1 


FW star 2 










(mA) 




26.00* 


490.514 


-2.050 


FMW 


73.63 


30.33 


26.01 


492.393 


-1.320 


FMW 


— 


— 


26.00 


492.477 


-2.220 


FMW 


107.96 


183.82 


24.00 


493.634 


-0.340 


MFW 


110.10 


82.99 


26.00* 


496.258 


-1.290 


FMW 


52.201 




28.00 


501.094 


-0.870 


FMW 


— 


— 


26.00 


504.422 


-2.150 


FMW 


169.70 


181.76 


26.00 


504.982 


-1.420 


FMW 


— 


— 


26.00* 


505.465 


-2.140 


FMW 


25.21 


25.79 


22.00* 


506.406 


-0.270 


MFW 


2.45 


2.50 


26.00* 


508.334 


-2.958 


FMW 


1.16 


11.24 


26.00 


509.078 


-0.400 


FMW 


— 


4.39 


28.00* 


509.441 


-1.080 


FMW 


23.30 


22.79 


26.00* 


512.735 


-3.307 


FMW 


72.89 


74.44 


26.00* 


514.174 


-2.150 


FMW 


73.17 


91.99 


26.00* 


514.373 


-3.790 


FMW 


27.32 


28.61 


22.00* 


514.547 


-0.574 


MFW 


— 


— 


26.00* 


515.191 


-3.322 


FMW 


83.65 


104.74 


26.00* 


515.906 


-0.820 


FMW 


74.80 


80.40 


26.00 


518.006 


-1.260 


FMW 


— 




26.00* 


518.791 


-1.260 


FMW 


56.57 


37.12 


26.00 


519.495 


-2.090 


FMW 


73.31 


113.68 


26.00* 


519.547 


0.018 


K94 


85.83 


21.61 
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26.00 
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K94 




26.94 
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FMW 






26.00 
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-0.840 


FMW 


114.71 


172.90 


26.01* 
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-4.250 


K88 


57.75 




26.00 
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-0.970 


FMW 


129.66 


135.07 
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-0.907 


MFW 


61.55 


19.62 


24.00 
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-1.400 


MFW 


214.13 


203.87 


27.00* 
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0.060 


FMW 
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26.00 


536.162 


-1.430 


FMW 


269.68 


225.63 


26.00* 


537.371 


-0.860 


FMW 


82.82 


97.37 


26.00 


537.958 


-1.480 


FMW 


102.03 


113.47 


26.00 


538.634 


-1.770 


FMW 




- 


26.00* 


538.948 


-0.410 


FMW 


20.50 


85.97 


28.00 


539.233 


-1.320 


FMW 


— 


4.94 


26.00 


612.025 


-5.950 


FMW 


112.24 


113.17 


26.00* 


615.938 


-1.970 


FMW 


26.72 


40.58 


26.00 


622.674 


-2.220 


FMW 


23.82 


6.31 


26.01* 


636.946 


-4.253 


K88 


— 


— 


26.00* 


639.254 


-4.030 


FMW 


64.58 


67.31 


21.01 


660.460 


-1.480 


MFW 


142.65 


127.39 


26.00* 


672.536 


-2.300 


FMW 


61.93 


8.42 


26.00* 


848.198 


-1.647 


K94 


39.32 


30.20 


14.00* 


850.222 


-1.260 


KP 


14.75 


1.58 


26.00 


851.510 


-2.073 





152.76 


132.81 


26.00 


858.225 


-2.133 





146.53 


141.05 


26.00* 


859.295 


-1.083 


K94 


83.14 


81.46 


14.00* 


859.596 


-1.040 


KP 


42.03 


38.56 


14.00* 


859.706 


-1.370 


KP 


60.52 


47.79 


26.00* 


859.882 


-1.088 





71.36 


59.76 


26.00* 


860.707 


-1.463 


K94 


48.49 


49.34 


26.00 


861.180 


-1.900 


FMW 


217.20 


226.04 


26.00 


862.160 


-2.321 





128.62 


135.64 


26.00* 


867.474 


-1.850 


FMW 


46.99 


25.25 


14.00 


868.635 


-1.200 


KP 




4.67 


26.00 


868.862 


-1.212 


FMW 


351.04 


340.18 


26.00* 


869.870 


-3.433 


K94 


86.81 


87.78 


26.00* 


869.945 


-0.380 





77.67 


69.31 


26.00* 


871.039 


-0.555 


K94 


87.34 


86.76 


14.00* 


872.801 


-0.610 


KP 


24.91 


22.03 


26.00* 


872.914 


-2.951 


K94 


91.38 


96.13 



Notes, f"' http://kurucz.harvard.edu/LINELISTS/LINES/gfall.ref 



this code predict temperature structures consistent with plane- 
parallel ATLAS9 and ATLAS 12 models as well as spherically 
symmetric PHOENIX (Hauschildt et al.|[T999) and MARCS 



( Gustafsson et al.|2008[) models. Furthermore, [Neilson & Lester] 



( 2008| l showed that spherical model atmospheres predict inten- 
sity distributions that fit interferometric observations of red giant 
stars ( |Wittkowski et al.||2004l |2006a|b| l with center-to-limb in- 
tensity profiles from model atmospheres, and they determined 
stellar parameters consistent with results using ATLAS9 and 
PHOENIX models. 

For the A3 analysis, models were computed assuming plane- 
parallel geometry and using the opacity distribution functions to 
minimize computing time. Derivation of stellar parameters for 
the two stars was done in the following manner. First, synthetic 
colours were computed from a grid of model atmospheres span- 
ning a range in Teff of 3000 to 8000 K, and in log g of to 3 
with solar metallicity. Comparing the synthetic colours to the 
given colours, values for T^ff and log g were estimated. Next, a 
new grid of stellar model atmospheres and synthetic spectra for a 
range of Tgf[ and log g about the preliminary estimates was com- 
puted such that Areff = ±400 K and Alogg - +1, while also 
varying the iron abundance. Line data were taken from Kurucz 



databas^^ For each synthetic spectrum, a x^-f^^ was computed 
and a new best-fit T^s, log g, and [Fe/H] was determined. Using 
these values, a new grid was computed, varying T^ff, log g, and 
the silicon abundance. A new value for Teff log g, and the silicon 
abundance was found and the process was repeated for oxygen, 
carbon and calcium. 

The best-fit stellar parameters for Star 3 and Star 4 result- 
ing from the A3 analysis can be found in Tables [14] and [15] 
The results for Star 4 are consistent with the parameters given 
in Table ]2] However, the results for Star 3 do not agree. This 
disagreement is due to the method which was employed. This 
method did not use any specific absorption lines to constrain the 
gravity or abundance. Instead, the parameters were determined 
using ablind;(f^-fit. Furthermore, possible degeneracies between 
log g and various abundances were ignored. For instance, a syn- 
thetic spectrum for a model with T^ff = 4200 K, log g = 2.25, 
[Fe/H]= 0.6 and C/0= 0.15 had a fit to the Star 3 spectrum with 
a value that was < 10% different than the value for the 
best-fit model. 



http://kurucz.harvard.edu 
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3.2.4. A4 

A4 (R. Peterson) analyzed the optical spectrum of Star 1, as part 
of an ongoing analysis of standard stars spanning a wide range 
of temperature, gravity, and metallicity. For these analyses, stel- 
lar parameters and abundances were derived by matching each 
stellar spectral observation to theoretical spectra calculated with 
an updated version of the |Kurucz| ( 



and the static, one-dimensional ste lar atmosphere models se- 
lected from the grid of Castelli & Kurucz (2003 ). A4 interpo- 



T993bl l SYNTHE program 



lated an appropriate model for each star, and used as input a 
list of molecular and atomic line transitions with species, wave- 
lengths, energy levels, q/-values, and damping constants. 

The Kurucz ^ijrj lists of atomic lines with known energy 
levels ("laboratory" lines) were modified by comparing calcu- 
lations to echelle spectra of standard stars. Moving from weak- 
lined to stronger-lined stars, first each spectrum was calculated, 
then the gf- values were adjusted individually for atomic lines 
and as a function of band and energy for molecular lines, and, fi- 
nally, a guess was made on the identifications of "missing" lines, 
those appearing in the spectra but not in the laboratory line list. 
This process was iterated until a match was achieved in each 
case. Peterson! ([2008 ) shows an example of the fits achieved in 
the near-UV for tumoff stars, from solar to extremely low metal- 
licities. Fig.|5]shows optical spectra for stronger-lined stars. 

Stellar parameters were derived from the spectra, and not 
from colours. The effective temperature Tetf was constrained 
by demanding that the same abundance emerge from low- and 
high-excitation lines of the same species (usually Fe I), and by 
fitting the Balmer line wings in stars of Tgff ~ 5000 K or hot- 
ter The gravity was inferred from the wings of other strong 
lines; comparing Fel and Fell abundances provided a check. 
Demanding no trend in abundance with line strength set u^i^-. The 
iron abundance and other elemental abundance ratios stemmed 
from matching relatively unblended weak lines. The resulting 
uncertainties are typically 0.1 - 0.2 dex in [X/Fe] for element X. 
Peterson et al. ( [2001 ) l provide details. 

A4 first compared Star 1 against the metal-poor K1.5III 
giant Arcturus (a Boo, HR 5340, HD 124897) and the super- 
metal-rich K2,III giant n Leo (HR 3905, HD 85503), two well- 
observed stars. This quickly established Star 1 to be cooler than 
any giant with solar metallicity higher than one-third solar that 
A4 has analyzed before. Moreover, throughout the red region the 
star exhibited a multitude of absorption lines not seen in either 
of the other two K giants. The majority proved to be TiO lines. 
Consequently, the TiO line list of Schwenke ( 1998) , downloaded 
from the Kurucz websitj'^was added. 

Several iterations were required to fit reasonably well the ob- 
served spectrum for Star 1. The first iterations refined Tetf, and 
simultaneously [Fe/H], then log g and y^ic. Ultimately the best fit 
was obtained for a model with the parameters given in Table 12 
The fit does deteriorate below 6000A, where "missing" lines re- 
main significant. 

With one exception, there was no need to alter the relative 
abundance of any element with respect to that of iron in Star 1 
from the values adopted for the Sun. For nitrogen in Star 1, how- 
ever, the abundance was lowered to [N/Fe] = -0. 1 dex to match 
the multitude of CN lines in the red. That oxygen in Star 1 is 
solar, [O/Fe] - 0, was confirmed from three separate diagnos- 
tics: the [O I] lines at 6300.3 A and 6363.8 A, the high-excitation 
OI triplet at 7771.9 A, 7774.2 A, and 7775.4 A, and a fit of the 

' ' http://kurucz.harvard.edu/LINELISTS/GFHYPER100/ 

http://kurucz.harvard.edu/molecules/TiO/tioschwenke.idasc-gz 



red TiO lines. For the Sun the older, higher oxygen abundance 
log(0/H) = -3.07 was adopted. 

For Arcturus and ji Leo, relative nitrogen abundances were 
increased by 0.1 dex. Relative abundances were also changed 
for other elements, notably sodium and aluminum, the light 
elements Ca, Mg, Si, and Ti, and elements beyond the iron 
peak, whose proportions are all known to vary among old stars 
dSneden et al.|2008| l. 

The degree to which the A4 calculations match the observed 
spectra of Star 1 and other standards is illustrated in Fig. |5] 
which shows a comparison in two wavelength regions of the ob- 
served spectrum to the calculated spectrum for four stars: the 
Sun, Arcturus, Star 1, and /i Leo. Both wavelength regions are 
relatively free of TiO absorption in Star 1, but its inclusion is 
nonetheless important to better define the continuum. Both have 
many CN features, which must be closely approximated to de- 
fine both continuum and blends. 

Although the match is not perfect, it is largely satisfactory 
in all four stars. Whenever a line is observed to be significantly 
too strong in one star, it is usually also significantly too strong 
in the others - indicating either a missing or wrong identifica- 
tion, or an erroneous gr/-value. The Call core mismatch is due 
to the chromospheric contribution to this line. The log g values 
are supported by the agreement between Ca 1 and Ca 11 lines and 
by Fell lines. Temperatures are supported by the wide range of 
lower excitation potentials spanned by lines of Si I, Ti I, and Fe I. 



3.2.5. A5 

The A5 analysis was performed by A. Goswami. The stellar at- 
mospheric parameters (Teff, \ogg, and [Fe/H]), for Star 1 were 
determined by an LTE analysis of the equivalent widths of 
atomic lines using a recent version of MOOG ( Sneden|[l973 i. 
The 59 cleanest lines of Fe I and two lines of Fe II within the 
three wavelength ranges recommended for Experiment 1 were 
used in the analysis (see Tab. 10 1. The Fe lines covered a range 



in excitation potential (1.0 - 5.0 eV) and equivalent widths (10 
- 165 mA). The excitation potentials and oscillator strengths of 
the lines were from various sources listed in the atomic spec- 
tral line database from CD-ROM 23 of R. L. Kuruc^U Model 
atmospheres were selected from the Kurucz grid of model at- 
mospheres available at the Kurucz web sita^ Kurucz models 



both with and without convective overshooting (see § 5.6 1 were 
employed, the latter from those that are computed with better 
opacities and abundances and labelled with the suffix "odfnew". 
It was found that the derived temperatures from both the options 
agree well within the error limits. 

The microturbulence Umic was estimated at a given effective 
temperature by demanding that there should be no dependence 
of the derived Fe I abundance upon the equivalent widths of Fe I 



lines (see Tables [12] and 13 1. The effective temperature was ob- 



tained adopting the derived value of the microturbulence param- 
eter by the method of excitation balance, forcing the slope of 
the abundances from the Fe I lines versus excitation potential to 
be near zero. Using the Fe I/Fe 11 ionisation equilibrium, the sur- 
face gravity of the star was obtained. The results are listed in 
Table [12 



http://www.cfa.harvard.edu/amp/ampdata/kurucz23/ 
http://cfaku5.cfa.harvard.edu/ 
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Fig. 5. A4 analysis: The plot compares observed spectra (heavy lines) to calculated spectra (light lines) for Star 1 and three standard 
stars. The two panels display two separate wavelength regions, as indicated in A at the bottom of each panel. The stars are displaced 
vertically for clarity. Each y-axis tick denotes 10% of the full height of the normalized spectrum. The star name appears in bold 
above the continuum of each plot, followed by the model stellar parameters adopted for the calculation: effective temperature Tf-g, 
surface gravity logg', logarithmic iron-to-hydrogen ratio [Fe/H] with respect to that of the Sun, and microturbulence y^ic. Above 
the upper boundary of each panel appear identifications of the strongest lines in the calculation for fi Leo: the three digits after the 
decimal point of the wavelength in A of the line is followed by a colon for a "missing" line, and then the species identification. 



3.3. Phoenix model atmospheres (P) 



A m odel fit using Version 15 of PHOENIX (Hauschildt & Baron] 
|1999j ) was provided by C.I. Short. A grid of spherical LTE atmo- 
spheric models was computed with a parameter range as given 
in Table |4] The formal numerical precision for the fitting in Tf-g 
was about +30 K. The mass was held fixed at 1 M© to determine 
radii for each value of log g. Convection was treated in the mix- 
ing length approximation with a mixing length parameter, I, of 
1 .0 pressure scale height. Atomic line data were taken from the 
ATLAS9 list (Kurucz 1993a). Most molecular line data are from 



Rothman et al. (2005) and Kurucz (1993a) except for the line 
lists of CO ( |Goorvitch|1994| i and Tid lSchwenke|1998l ). 

For the purpose of Experiment 1 the approach of investigat- 
ing the global fit to two wavelength regions ("visible" and "near- 
IR", see Fig.[T]i was adopted in the P analysis. The radial veloc- 
ity correction was determined by fitting the Doppler shift, AA, 
to relatively unblended weak spectral lines at the blue and red 
ends of the observed spectra. Synthetic spectra were computed 
with a spectral resolution, R, of 300 000. Assuming that giants 
are slowly enough rotating that vsini would be minor compared 
to the effect of microturbulence, no rotational broadening was 
taken into account. An initial rectification was performed by 
dividing the synthetic spectra by unblanketed spectra for each 
model; the spectra were then re-rectified to the observed spectra 
piece-wise with a single-point calibration in each of the two fit- 
ting regions. The visible region is affected by significant telluric 
contamination. In principle, every feature in the entire spectrum 
contains information about the stellar parameters, provided they 
can be modelled correctly, and this experiment can be seen as a 



test of how well a model "blindly" fitted to the broad spectrum 
with equal weight on all spectral lines would recover the stellar 
parameters of a known star. 

In the near-IR region the strongest diagnostics are the Ca II 
IR triplet lines. However, it was found that PHOENIX con- 
sistently produces profiles for these lines that are too bright 
throughout the damping wings for any realistic stellar parame- 
ters, even with line damping parameters tuned to match the solar 
line profiles. A comparison of unblanketed spectra among the 
workshop participants revealed that PHOENIX predicts a larger 
"background" continuum flux throughout the near-IR band than 
the other codes, and a possibility to be investigated that is consis- 
tent with both of these results is that PHOENIX under-estimates 
the continuous opacity in this region. 

Figure|6]shows the relative flux residuals, (/i, o s)//i, o, 
where f^,o and /^,s are observed and synthetic flux, respec- 
tively, for Star 1 and models of [M/H]= 0.0. Figures |7] and |8] 
show the values of P^r interpolated A point for Star 1 and 
Star 2, respectively, and all models. 



Star 1 For the visible fitting region, the values showed dis- 
tinct global minima as functions of Teff (minimum ;t'^(A/l)"' » 
0.025), but picked out best-fit values of log g and [M/H] less 
clearly. Models of [M/H]=0.0 yielded best-fit T^^ values of 3940 
and 3875 (+30 K, formally) for log g values of 2.0 and 1.5, re- 
spectively. For [M/H] =-0.5 these \ogg values yielded best fit 
Teff values of 3800 K. Models of loggr=1.0 provided a signif- 
icantly worse fit. There was marginal evidence that models of 
[M/H]= -0.5 might provide a better fit than those of [M/H] = 
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Table 10. Lines analysed by A5 (Sect. |3.2.5 1. For each line, the 
wavelength, the gf value, the lower level excitation potential 
E\ovi, the equivalent width (EW), and the source of the Une data 
according to |Kurucz & Bell|fl995| )'' are given. 
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-1.665 
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Notes. http://kurucz.harvard.edu/LINELISTS/LINES/gfall.ref 
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Fig. 6. P analysis: Relative flux difference, {f^^o - /i,s)//i,o, 
smoothed to a resolution corresponding to A/l - 50 A, for a 
subset of the grid with [M/H]=0.0 around the best-fit parameters 
for Star 1 are shown in gray-scale. The best-fit model (Teff/log 
g/[M/H] = 3940/2.0/0.0) is overplotted with a thicker line style. 



0.0. For the near-IR region the values showed a much more 
shallow minimum around 3800 K (minimum x\^'^T^ ~ 0.01) 
for [M/H]= 0.0 and log g values of 2.0 and 1.5, and did not show 
a minimum at all for [M/H]= -0.5. 

Star 2 The fits yielded a similar result as that for Star 1, but 
with all \hQ x^i^^y^ minima shifted downward in T^g by about 
100 K for the visible band, and 50 K for the near-IR band. For 
the visible band the best fit value for log g was either 1.5 or 2.0, 
and the best fit values for T^s were 3800 K ([M/H]= 0.0) or 
3700 to 3750 K ([M/H]= -0.5). Again, only models of [M/H] = 
0.0 showed a minimum for the IR band, with a best fit Tgff value 
being indicated by a very shallow minimum of ~3750 K. Models 
of [M/H] of 0.0 and -0.5 provided the same quality of fit, so the 
nominal best fit value was found to be -0.25 + 0.25. 

In conclusion, the near IR band needs attention in the P mod- 
elling to cultivate it as a more effective diagnostic of stellar pa- 
rameters for M giants. The lack of distinct minima in the fitting 
statistic may be caused by the poor quality of the fit to the broad 
Ca n IR triplet lines, as discussed above. Therefore, the best-fit 



parameters for the P analysis given in Tables 12 and 13 are based 
entirely on the fit to the visible region. 

3.4. CODEX model atmospheres (C) 

The authors contributing the C analysis were M. Ireland and 
M. Scholz. The Cool Opacity-sampling Dynamic Extended 



(CODEX) models as published inllreland et al. (2008 1 are based 



on a combination of three separate codes - a grey dynami- 
cal atmosphere code, a temperature iteration code and a spec- 
trum computation code. In addition, the equation of state is pre- 
tabulated based on the 1970s code from T. Tsuji, including up- 
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Fig. 7. P analysis: Values of;(f^(A/l)"' for Star 1. Solid lines: Fit 
to visible region; Dotted lines: Fit to near-lR region. No point 
symbol: [M/H]= 0.0; Asterisks: [M/H]= -0.5. Dark to light gray 
lines: log g- 2.0 to 1.0. 
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Fig. 8. P analysis: Same as Fig.jrjexcept for Star 2. 



dates from'Sharp & Huebner'(T990| and|Jeong et al.'('200T). The 
CODEX models are designed to predict the observables in the at- 
mospheres of Mira variables, and the main algorithms involved 
are described in |Schmid-Burgk| ( |197 5 ) and Schmid-Burgk &| 
Scholz (1984 1. For the static models computed for this paper, 
the dynamical code was not used and the pressure stratification 
was iterated simultaneously with the temperature. 



An important property of these models relevant to the test 
stars in this paper is that only atomic lines of neutral atoms 
were included, so that the code was not applicable and has not 
been used for effective temperatures above 4000 K. Furthermore, 
there were numerical instabilities that are not fully understood, 
which means that the spherical code could only be used in sig- 
nificantly extended configurations. In practice, this implied a 
restriction to loggf values less than approximately -0.5 to 0. 
Within these constraints, the model structures produced by the 
CODEX models agree with MARCS structures within 50 K. The 
H20 lines come from Partridge & Schwenke ( 1997[ l and those 
of TiO from Schw enke] ( |199 8V Other diatomic molecular (CO, 
OH,CN, SiO, MgH) and neutral metal atomic lines come from 
the input to the ATLAS 12 models ( |Kurucz||1994[ l. VO was not 
taken into account here. 

The full grid of models used to fit the stars in Experiment 1 
had the parameters given in Table |4] and a fixed extension 
of 100 times solar. The atmospheric extension in this context 
was defined as (reif/5770)(/?/Ro)(M/Mo)-'. A spectral fit was 
obtained at fixed gravity by minimising the mean square dif- 
ference in between the model and the sample spectrum after 
interpolating bi-linearly in metallicity and effective tempera- 
ture. Additional fit parameters were macroturbulence (or almost 
equivalently spectrograph resolution/stellar rotation velocity), 
radial velocity, continuum flux and continuum slope. A single 
continuous spectral window was used for the analysis of Star 2 
(see Fig.[T]). 

The only stars for which reasonable spectral fits could be ob- 
tained were Star 2 and Star 4, as they were the only stars with 
parameters that were covered by the model grid. Uncertainties in 
the fitted parameters were not estimated in the C analysis, but fits 
were noticeably poorer by eye at parameters different by 0.2 dex 
in metallicity or 50 K in effective temperature, and so these can 
be taken as uncertainties. Gravity was not meaningfully con- 
strained, but the log g value with the marginally better fit was 
chosen from the two available. The best-fit parameters are given 
in Tables 13 and 15 For Star 4, the broad band colours computed 
from the models are V - 7=1.89, J - K^Q.9A and V - K^3.1\. 



3.5. TsujiOS model atmospheres (T) 

The T analysis was provided by T. Tsuji. Modelling was re- 
stricted to Star 4, as the recent focus of the contributing author 
has been on the analysis of near infrared molecular lines ( Tsuji 
|2008[ ). With the given photometric data (see Table |2]i it was not 
easy to infer fundamental parameters of this star According to 
Perrin et ar] ([T998 ) J-K - 1.21 corresponds to an M7 giant (see 
their Table 3), which would have a typical Teff of 3087+94 K 
(see their Table 5). Tgff - 3100K was, thus, chosen for Star 4. 



In the absence of any information on mass, radius, etc. it was 
simply assumed that M - 2Mq and R - 2OOR0. These model 
parameters were adopted, together with a microturbulence i;n,ic 
— 3km s ' , from a grid of model photospheres for red giant stars 
presented in |Tsuji| ( |2008j ). The abundance pattern used for the 
modelling is described in Table 4, case (b), in Tsuji ( 2008| l for 



CNO and in Table 1 in Tsuji (2002i for other elements. Line 



data for CO 3-0 lines are listed in Table [TT| Line positions and 
gf values are referenced in (Tsuji 2008 ). In particular, log Ac/Ah 
= -3.92, log An/Ah = -3.55, log Aq/Ah = -3.31, and [Fe/H] = 
0.0, were chosen. 

For the abundance analysis, 14 CO lines relatively free from 
blending in the given spectrum of Star 4 were measured. The 
resulting equivalent width (EW) data are given in Table 11 A 



measurement of OH and CN lines was also attempted but it 
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Table 11. Measured CO 3-0 lines in Star 4 from the T analysis 



(Sect. 3.5 1. For each line, the wavelength A (air), the gf value, 
the lower level excitation potential Ziiow, the equivalent width 
over wavenumber log W/v, and the identification are given. 



A 


log(g/) 


^low 


log W/v 


Identification 








(nm) 




[cm-'] 








o 






1560.082 


-7.797 


886.903 


-4.696 


'^C 


^0 3-0 R 21 






1560.447 


-7.823 


806.383 


-4.632 


'^C 


^0 3-0 R 20 


o 






1560.839 


-7.850 


729.678 


-4.691 




^0 3-0 R 19 


< 






1561.257 


-7.877 


656.789 


-4.679 


'^c 


^0 3-0 R 18 






1561.533 


-7.249 


5055.604 


-4.732 




^0 3-0 R 51 








1561.702 


-7.907 


587.721 


-4.621 


'^c 


^0 3-0 R 17 








1562.002 


-7.234 


5252.117 


-4.784 


'2c 


^0 3-0 R 52 








1562.499 


-7.220 


5452.275 


-4.744 


'^c 


^0 3-0 R 53 








1563.748 


-8.038 


349.698 


-4.645 




^0 3-0 R 13 








1564.157 


-7.179 


6074.532 


-4.791 


'^c 


^0 3-0 R 56 








1564.326 


-8.076 


299.766 


-4.695 




^0 3-0 R 12 




LT) 




1564.930 


-8.116 


253.667 


-4.707 


'^c 


^0 3-0 R 1 1 


< 


d 




1565.561 


-8.159 


211.404 


-4.729 




^0 3-0 R 10 


Alog 


o 




1566.763 


-7.125 


6954.739 


-4.902 


'^c 


^0 3-0 R 60 


















d 





was diflicult to find a suflicient number of lines free from dis- 
turbance. Line-by-line (LL) analysis as described in detail in 
|Tsuj i (1986, 2008 ) was applied to measure the selected CO lines. 
Results are shown in Fig. |9^, in which logarithmic abundance 
corrections to the assumed carbon abundance for different lines 
are given for ii^ic = 3, 4, and 5 kms"'. For the micro turbu- 



lence value given in Table 15 i, the abundance corrections were 
found to be independent of EWs, and the mean logarithmic abun- 
dance correction was found to be 4-0.36 + 0.13. This resulted in 
logAc/AH = -3.92 + 0.36 = -3.56 + 0.13. Fig.^ shows thatthe 
logarithmic abundance corrections are on average zero for these 
parameters. 

For the resulting microturbulence and carbon abundance, 
a synthetic spectrum for the region of the CO 3-0 band was 
computed, and the result is compared with the observed one in 
Fig. [TO^ . The CO lines used in the LL analysis are indicated by 
the rotational identifications for the CO 3-0 lines. Only these 
CO lines should be considered for estimating the quality of the 
fit, but fits are not so good even for these lines. To improve the 
fits in Fig. 10?, some additional broadening may be required, 
and several values of macroturbulent velocities were tried. It was 
found that the fits are somewhat improved for the selected CO 
lines with an assumed macroturbulent velocity of 3.0 km s-1 
(which could have been determined empirically if a sufficient 
number of weak lines could be measured, as in Fig. 3 in 
([1986)). 



Tsuji 



Because the stellar parameters (Teff , M, /? etc.) cannot be pre- 
defined well and because the number of lines is too small for 
the LL analysis, only the basic procedure to analyze the stellar 
spectrum has been shown and a final converged solution was not 
pursued. However, we may conclude that the abundance deter- 
mination is possible even for very cool stars in this way, if dozens 
of weak lines could be measured. 

For the T analysis, the classical theory of line formation as- 
suming the presence of microturbulence was applied. However, 
the applicability of the classical model of line formation may 
be limited to relatively weak lines with log W/v < -4.75 (Tsuji 
2008[ l. Although the sample of CO lines of Star 4 analyzed 



here includes some lines slightly stronger than this limit, the 
modelling was quite successful because most of the lines used 
are weak and below or near the limit mentioned (|Tsuji 1991 



cf. a similar previous study of the CO second overtone 
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Fig. 9. T analysis: a) Line-by-line analysis of CO lines in Star 4, 
i.e. logarithmic abundance corrections to the assumed carbon 
abundance derived from the observed log W/v values are plot- 
ted against the log W/v values for Umic = 3, 4, and 5 kms^^ b) 
Confirmation of the null abundance corrections for logAc/Au = 
-3.56 and Um\c = 3.93 km s"' . 



The strongest lines with log W/v > -4.5 clearly deserve an 
extra contributions from outside the photosphere ( |Tsuji|[T988] l. 
Thus, if a large number of the intermediate strength lines with 
-4.75 < log W/v < -4.5 are included as in a previous study of 
the CO first overtone lines (Tsuji 1986]), the results cannot be 
correct as was recognized 20 years later ( [Tsuji 2008 1. Such a 
problem will appear not only in very cool stars but also in earlier 
type stars such as K giants ( [Tsuji|2009[ ). 



ines). 



4. Results 

As described above, two experiments were conducted. 
Experiment 1 had the aim to find the best fitting model to the 
four sample spectra provided, using the preferred method of each 
participating researcher or team. It was to be completed before 
the beginning of the workshop. A time limit of several weeks 
was set for the groups to deliver their results. This limit was in- 
troduced not only for practical reasons but also to simulate the 
typical amount of time used for parameter determination from a 
stellar spectrum in the current age of multi-object spectrographs 
and large surveys. Naturally, this time limit set some constraints 
on the level of elaboration of the resulting fits. Experiment 2 re- 
sulted from the discussions at the workshop and was aimed at 
a direct comparison of synthetic spectra from different model 
codes for a given set of stellar parameters. 

4.1. Experiment 1 - stellar parameter determination 

We start with an overview of the main aspects of the differ- 
ent modelling approaches (for details see S ections [3 . 1 1 to [33]l . 



18 



T. Lebzelter et al.: Comparative Modelling of the Spectra of Cool Giants 



a) star4/CO/K^^^ = 0,0 
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Fig. 10. T analysis: a) Observed spectrum (dots, relative flux vs. 
wave number) and model spectrum (solid line) with a macrotur- 
bulence of zero, for the CO lines between 1.5601 and 1.5625 yum 
of Star 4. b) The same as a) but with a macroturbulence of 
3.0 kms'. 



About half of the groups (six of fourteen) based their analy- 
sis on MARCS atmospheric models calculated in spherically- 
symmetric geometry. Five groups used the plane-parallel 
ATLAS suite of codes ( Kurucz 1993a b , Lester & Neilson 2008f 
The three software packages PHOENIX ( Hauschildt & Baron] 
T999] l, CODEX preland et al.p OOS ) and TsujiOS ( |Tsuji||2008| l 
were employed by one group each (all of them spherical). The 
first character of the group identifier indicates this aspect. The 
stellar radius, a key parameter, is defined as the layer where 
the Rosseland optical depth is equal to one in the case of the 
MARCS models, the model of Tsuji, and the CODEX models. 
In the case of the PHOENIX models, it is assumed that the stellar 
mass M is one solar mass, and the radius R is computed accord- 
ingly from the specified surface gravity (R - (GM/gy^^). For 
the plane-parallel ATLAS models the radius is undefined. 

To improve the fit of the line widths and to take into account 
various additional line broadening mechanisms like stellar rota- 
tion a macroturbulence has been added by most of the groups, 
i.e. the synthetic spectra have been convoled with a Gaussian 
profile. 

The set-up did not require all participants to complete the 
experiment for all four spectra. Indeed, due to methodological 
constraints, some participants restricted their attempts to either 
the visual (M2, M3, M4, M6, A2, A4, A5, P) or the near in- 
frared spectra (M5, A3, T). Naturally, the team that prepared the 
artificial observations for Star 3 and Star 4 (Ml) did not partic- 
ipate in the fitting of these spectra. Two groups analysed both 
wavelength domains - Al {Star 1 and Star 3) and C {Star 2 and 
Star 4). In summary. Star 1 was analysed by eleven groups, of 
which five also analysed Star 2. Star 3 and/or Star 4 were anal- 
ysed by five groups. 



4.2. a Tau CStar ^)anda Cet fStar 2) 

The best fits to these two spectra as provided by the participating 
teams are summarized in Tables [T2] and [T3] For a Tau we have 



the more extensive set of model fits. The results are presented in 
Fig. 1 1 for the three possible combinations of the three main pa- 
rameters derived: T^ff, log g, and [Fe/H]. The (unweighted) mean 
values of all results are: T^s = 3980+250 K, log g = 1 .3±0.3, and 
[Fe/H] = -0.2+0.2 dex, close to the literature values discussed 
in Section 2.2 and Table [T] No systematic differences between 
the model families (MARCS, ATLAS, and others) are apparent. 
Within the error bars given by the different groups, the derived 
Teff values of almost all groups agree with each other The only 
exception is the A2 analysis, which deviates by 800 K from the 
mean value. Discarding the A2 value, the mean Tgff of a Tau 
as a result of this experiment becomes 3910±80 K. Note that in 
the case of our experiment, contrary to a series of repeated mea- 
surements, the mean value does not necessarily result in a more 
precise value than the individual determinations. 

We have no obvious explanation for the deviating result of 
A2, but the analysis approach differs in various aspects from the 
bulk of the other analyses. First, A2 uses a large model atmo- 
sphere grid including T^ff values up to 6000 K, which makes 
it possible to probe higher temperatures than other analyses. 
However, the same can be said for the M3a and M3b analyses, 
which do not find a high-Tetf solution using a grid up to 8000 K. 
Second, A2 did not include molecular lines in the analysis, but 
the same is true for M6, Al, and A5. A final aspect is that A2 
did an equivalent width analysis, and not a synthetic spectrum 
fit, which seems problematic when dealing with crowded spec- 
tra. The A5 analysis also used equivalent widths, and obtained a 
Teff close to the mean value. However, there are important differ- 
ences between the two analyses, which might explain the differ- 
ent results. A2 measured equivalent widths for selected lines of 
several different species from the observed and computed spec- 
tra, while A5 computed equivalent widths line-by-line, adjusting 
the abundances for lines of Fe I and Fe II only. Also, A2 derived 
all parameters simultaneously from a fit to all lines, while A5 
used the classical approach of excitation and ionization equilib- 
rium for Fe lines. Note that there are only three lines in common 
in the line lists of A2 and A5. 

Concerning gravity, the results for a Tau group around two 
values, slighdy lower (Ml, M3b, M6, A4, A5) and sHghtly 
higher (M2, M3a, Al, A2, P) than the published values given 
in Table [T] with the two M4 results falling in both groups. The 
majority of the modelers found a Tau to exhibit a lower metal- 
licity than the Sun, the lowest values being suggested by Al. 
For the cases where [Ca/Fe] values were also determined, most 
groups find a mild underabundance of Ca. 

Figure [12] shows the best-fit T^ff and log g for Star 1 {a Tau) 
from all participants except A2 converted to luminosity, using 
three different mass values, together with the direct T^s and log g 
values from Table]!] We also show evolutionary tracks for metal- 
poor models by the Padova group for several different mass val- 
ues. It is obvious that some of the derived pairs of parameters are 
inconsistent with the chosen stellar evolution models. 

Fewer model fits are available for a Cet, and thus we have 
omitted a graphical presentation. The (unweighted) mean val- 
ues and standard deviations of all results are: T^.^ - 3830+220 K, 
log g - 1.0+0.8, and [Fe/H] = -0.2±0.2 dex, which is close to the 
published values discussed in Section 2.3 As for a Tau we find 
rather good agreement in the derived temperature, with the ex- 
ception of A2 (deviating by 500 K). Discarding the A2 value, 
the mean Teff of a Cet as a result of this experiment becomes 
3760±70 K. Concerning the surface gravity, the scatter is some- 
what larger than for a Tau, and the majority of the values are 
larger than the published value. For [Fe/H] the model fits to the 
spectrum of a Cet suggest metallicities from solar to subsolar In 
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Fig. 11. Results of the various model fits in three different two-parameter planes for Star 1 (a Tau). The red star marks the values of 
the direct parameters given in Table [T] 
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Fig. 12. Best-fit T^s and log g for Star 1 (a Tau) converted to 
luminosity, using three different mass values (triangles up - 
0.6 Mq, squares - 1.1 Mq, triangles down - 1.6 Mq), and direct 
Teff and log g values from Table [T] (black cross with error bars). 
Solid lines are evolutionary tracks for Z=0.008 and Y-Q.26 by 
the Padova group (Bertelli et al.|2008| [20091 ), labelled to the left 
with mass in solar masses. 



general, one would expect a correlation between T^f^ and [Fe/H] 
for the differing results, since line strength increases towards 
lower temperatures and higher [Fe/H]. This can be seen to some 
extent in Fig. [TT] and the a Cet results show a similar trend. On 
the other hand, there is no obvious correlation between log g and 
Teff, or log g and [Fe/H]. 



4.3. Star 3 and Star 4 

The fitting results for the two artificial stars are summarized in 
Tables [14] and [15] In this part of the experiment we asked the 
participants to derive also a C/O ratio. 

The fitting of spectra in the optical range is still more widely 
used than the application of the same technique to spectra in 
the H and K band. This may be due to the longer and wider 
availability of optical high-resolution spectrographs and the bet- 
ter knowledge of the resident features. It is therefore not surpris- 
ing that also the number of teams providing a fit to one of the 
near-infrared spectra is smaller than for the case of the optical 
ones. Furthermore, the lower temperature, especially of Star 4, 
requires an extensive inclusion of molecular line data, which is 
not needed for the usual work of some of our participating teams 
and, therefore, not included in their codes. 

In comparing the fit results with each other we have to re- 
member that the fitted star' was a MARCS model. Therefore, 
one could expect that systematic differences between the un- 
derlying codes of the model fits should show up more clearly. 
For Star 3 the result is a bit surprising, because the best fits of 
all three groups gave a solar or supersolar metallicity, while the 
Star 3 model is clearly subsolar. Also, the derived C/O ratio dif- 
fers significantly from the model input value. Temperature val- 
ues, on the other hand, nicely agree with each other and are in 
reasonable agreement with the Star 3 value. For Star 4, all model 
fits are close to the real values or at least point into the correct 
direction (in particular to a higher C/O ratio). 

A conclusion may be that the spectra of cool stars can be un- 
derstood reasonably well if the attention is confined to relatively 
weak lines. However, we are still far from understanding com- 
plete spectra of cool stars including the intermediate strength 
lines and the strong lines, and the status of the modelling of these 
strong lines is still not at a satisfying stage. It appears to be no 
problem in analyzing the CO lines in Star 4, but this may be 
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Table 12. Experiment 1 : parameters determined for a Tau (Star 1) from the best fits to the optical spectra by the various participating 
teams. 



starting parameters from photometry best fit parameters from spectroscopy additional quantities 



Code 




logg 


Lre/HJ 




logg^ 


rr7„ /TTi 


LCa/t<eJ 


KV 










{g (cm s"^)) 






(g (cm s"-)) 






(Km s ) 


(Km s ) 


(Km s ) 


Ml 






-0.3" 


3800±100 


0.75+0.25 


-0.3+0.25 




54 




y 


M2 


3850/4000'' 






3850 


1.25 


0.0 


0.0+0.1 


54.5 


1-2' 


4.5 


M3a 








3994±100 


1.54+0.20 


-0.15+0.10 




53.9+0.9 




9 












U.y4+U.l J 


— U.JJ+U.IU 


— U.Z+U. i 


jj.y+u.y 


of 


Q 


M4 








4000'' 


1.5'' 


0.0'' 




55 


2f 












3900'' 


1.0'' 


-0.25'' 




55 


If 




M6 


3750-4250'' 


1.0-1.5'' 




4000+200 


1.0+1.0 


0.0+0.5 


-0.3+0.1 


54.3+0.5 


If 




Al 


3900 


0.0 


1.0 


3950+100 


1.5 


-0.67+0.08 


-0.23+0.17 




2 


( 


A2 








4788±92 


1.5+0.1 


-0.4+0.1 




54.1+0.4 


If 




A4 








3900 


1.0 


0.0 




54.3 


1.3 


1.5 


A5 








4000+200 


1.0+0.25 


-0.29 






1.6+0.5 




P 








3910+45'' 


1.75+0.25'' 


0.0'' 




54.71 


3'- 













3800+30'' 


1.75+0.25'' 


-0.5'' 




54.71 


3" 






Notes. Photometric colours used to resolve a Teff [Fe/H] degeneracy after spectral fitting. reff(V - K)ITcff{J - K). [<i/Fe]. '''' Degeneracy 
in Teff and [Fe/H]. *''' Ranges. *'^' Fixed value(s). Not taken into account. Microturbulence increased from 2.0 to 4.0 kms"' as logg decreased 
from 2.5 to 1.0. Not determined. 

Methods for photometric parameter determination: Ml, M2: MARCS models. Theoretical colour {V - K and J - K) - T^g calibrations ( |Bessell| 
let al.|1998|Fan Eck et al.|20TT) . M6: |Houdashelt et al.|pOOOt . Al: |Kucinskas et al.|p005| l. 



Table 13. Experiment 1: parameters determined for a Cet from the best fits to the optical spectra by the various participating teams. 





starting 


parameters from photometry 


best fit parameters from spectroscopy 


additional quantities 


Code 


7'eff 


log 9 


[Fe/H] 




log 9^ 
(g (cm s"-)) 


[Fe/H] [Ca/Fe] 


RV 






(K) 


ig (cms-^)) 




(K) 




(kms-') 


(kms"') (kms"') 


Ml 






-0.5" 


3675+50 


0.5+0.25 


-0.5+0.2 




2.5'^ 3'' 


M2 


<3700 






3700 


1.0 


0.0 


-25.9 


1-2'^ 5.0 


M3a 








3867+100 


1.15+0.20 


-0.22+0.10 


-27.6+1.3 




M3b 








3718+80 


1.31+0.15 


-0.02+0.10 -0.32+0.10* 


-27.6+1.3 


2" '' 


Al 


3500 


1.0 


-2.0 












A2 








4310+136 


1.5+0.1 


-0.6+0.1 


-26.3+0.4 


2' 


P 








3750+30 


1.75+0.25 


-0.25+0.25 


-26.22 


y 


C 








3820+50 


-0.5 


0.0+0.2 




2.8'' 2.3 



Notes. *"' Photometric colours used to constrain metallicity after spectral fitting. [a/Fe]. Fixed value(s). '''' Not taken into account. 
*'"' Microturbulence increased from 2.0 to 4.0 km s"' as log g decreased from 2.5 to 1 .0. 

Methods for photometric parameter determination: Ml, M2: MARCS models. Theoretical colour (V - K and J - K) ~ T^ff calibrations ( |Bessell| 
let al.|1998|rvan Eck et al.pOTT) . Al: |Kucinskas et al.|pU05> . 



simply because the provided spectrum of Star 4 does not include 
strong lines, and we should not be fully satisfied with the present 
result. 



4.4. Experiment 2 - synthetic spectra comparison 

For Experiment 2 the participants were asked to calculate a syn- 



thetic spectrum using the stellar parameters listed in Sect. 2.4 



Spectra should be provided for the wavelength ranges 4900- 
5400 A, 6100-6800 A, and 8400-8900 A, respectively. Eight 
high-resolution spectra (R-300 000 to 500 000) were produced 
for the comparison, namely P, Al, A2, A3, Ml, M2, 1V13, and 
M4. Note that M3 is at that stage (no parameter determination) 
very similar to M2 since the spectra are based on the same line- 
lists and model atmospheres but due to interpolation and the 
selected parameters the spectra are slightly different. A3 pro- 
vided both a plane-parallel and a spherical solution, which are 
very similar. The plane-parallel version is used in the following 
as it had the higher spectral resolution. Al and A2 both used 
ATLAS9-I-SYNTHE, but under two different operating systems. 
They differ only at very few wavelength points, so that we show 



only Al in the figures below. An overview of the modelling de- 
tails can be found in Table [T6l 

To compare these synthetic spectra we first made a linear 
interpolation of the normalized flux values to a common wave- 
length grid with a step size of 0.01 A. Then, for each wave- 
length point the mean flux of all six spectra (P, Al, A3, Ml, 
M2, M4) was calculated, and the root-mean-square (rms) differ- 
ence of all normalized fluxes from the mean normalized flux was 
determined. Next, the rms differences were averaged over 2 A- 
wide wavelength intervals. In addition, the difference of each 
individual normalized spectrum from the mean normalized flux 
was calculated. 



Figure 13 shows the rms flux differences. We note that a large 



rms value can be caused by one or more of the following: a large 
difference for an individual line, small differences for several 
lines, or a small continuum difference. To explore the features 
seen in Figure 13 we looked at the individual spectra in more 
detail. In Fig. 14 we show the differences of the six spectra from 



the mean, for selected 10 A-wide wavelength intervals. Starting 
with Region 1 in Fig. [13] we found that MgH is the dominant 
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Table 14. Experiment 1: parameters determined for the artificial spectrum of Star 3 from the best fits of the various participating 
teams. 





starting 


parameters from 


photometry 




best fit parameters from spectroscopy 






Code 




logg 


[Fe/H] 




logg^ 


[Fe/H] 


C/O 








(K) 


(g (cm s"^)) 




(K) 


(g (cms"-)) 






(kms-') 


(kms"') 


M2 


4225 


















M5 


4250 


2.0 


0.0 


4250 


2.0 


0.0 


0.46 


1.7° 


4.5 


Al 


4300 


2.0 


-2.5 


4350±100 


2.0 


0.32±0.04 


0.63* 


2 


C 


A3 


4200 


3.0 


0.0" 


4400 


3.5 


0.6 


0.63 


2.0" 





Notes. Fixed values. <*' [C/H]=0.31, [N/H]=0.02, [O/H]=0.21. Not determined. 

Methods for photometric parameter determination: M2: MARCS models. Theoretical colour (V — K and J - K) - T^ff calibrations ( [Bessell et al.| 
[T998i|van Eck et al.|2011> . M5: |Worthey & Lee| ( |20i 1, 1. Al: ,Kucinskas et aL| ( |2005) . 

Table 15. Experiment 1: parameters determined for the artificial spectrum of Star 4 from the best fits of the various participating 
teams. 



Code 


starting parameters from photometry 
Teff log g [Fe/H] 
(K) (6f(cms-2)) 


(K) 


best fit parameters from spectroscopy 
log g [Fe/H] C/O 
{g (cm s-2)) 


(kms-') 


f'mac 

(kms-') 


M2 


>3300 K/3500 K" 


















M5 


3250 


0.0 


0.0 


3500 


0.0 


0.0 


0.69* 


3.0'- 


5.0 


Al 


3800 


1.0 


-1.5 














A3 


3000 


1.0 


0.0' 


3200 


0.25 





0.79 


2.0'' 





C 








3244±50 


0.0 


-0.1+0.2 


0.48 


2.8'' 


1.2 


T 


3087±94 






3100 


0.14 


0.0 


0.56 


3.9±0.4 


3.0 



Notes. Teff (V - K)IT^e(J - K). <*' C/O not consistent with model structure. '''' Fixed values. 

Methods for photometric parameter determination: M2: MARCS models. Theoretical colour {V — K and J - K) - T^g calibrations ( [Bessell et al.| 
[19981 [van Eck etal.|2011> . M5: |Worthey "& Lee (20lT). Al: |Kucinskas et aL| ( |2005) . T: |Perrin et af]fT998| l. 

Table 16. Overview of model inputs used for Experiment 2. The listed codes were also used for Experiment 1. 



Label 


Codes 


Spectral 
resolution 


Comments 




Ml 


MARCS + comaOS 


300000 


spherical 




M2 


MARCS 

+ TURBOSPECTRUM 


500000 


mass = 1 Mq; spherical 




M3 


MARCS 

+ TURBOSPECTRUM 


500000 


mass = 1 Mo, log9=1.5, [Fe/H]=-0.25 
spherical 




M4 


MARCS + BSYN 


300000 at 6000 A 


mass = 2 Mq, y,nic=2kms-'; spherical 
spectrum convolved to a resolution of 80 000 




Al 


ATLAS9 + SYNTHE 


500000 


model atmospheres interpolated within Kurucz 


grid 


A2 


ATLAS9 + SYNTHE 


300000 


model atmospheres and spectra calculated 






Linux vers. (Sbordone) 




for [Fe/H]=0 and -0.5, interpolated to [Fe/H] = 


-0.2 


A3 


ATLAS 12 


500000 


plane-parallel 




P 


Phoenix V. 15 


300000 







molecular absorber in the region from 5000 to 5210 A, where 
we see high rms differences. Comparison with a spectrum cal- 
culated with MgH lines only demonstrated that the differences 
in the spectra produced by different investigators are mainly due 
to differences in the MgH line data that they used. The broad 
rms peaks around 6220 A in Region 2a and around 6740 A in 
Region 2b coincide with regions where TiO absorption domi- 
nates, while the smaller peak around 6510 A in Region 2b coin- 
cides with dominant CN absorption. 

The three sharp peaks in Region 2a close to 6350 A are due 
to three atomic lines. Two lines at 6343 and 6362 A are present 
in four out of the six spectra, and are much broader in the A3 
model than in the others. These lines are not present in the ob- 
served spectrum of a Tau. One line at 6354 A is present in all 
models, but with varying line depths (0.3, 0.4, 0.6, 0.8, 0.9) - a 
comparison with the observed spectrum of a Tau points towards 



intermediate values. A careful check of the line lists revealed 
that these three lines can likely be associated with wrong data 
for some Cai transitions. A3 included theoretical iron line lists 
from Kurucz, which are based on predicted energy levels. These 
lines appear to cause the strong and broad features seen in A3's 
synthetic spectra at the mentioned wavelengths. A minor sharp 
peak can be seen in Region 2b close to the Ha line, which was 
not included in the line list of the M2/M3 model (see Fig. 14 i). 
We note that none of the models reproduce the Ha line in the 
observed spectrum of a Tau. The observed line depth is more 
than twice as large as any of the calculated ones. In Region 3 in 
Fig. 



13 



the most obvious features are centered on the IR Ca II 
triplet lines and are due to different modelling of the broadening 
of these lines. A3 and M2/M3 show the broadest lines, A1/A2, 
Ml, and M4 are narrower, and P is narrowest. An example is 
shown in Fig. 14 ;. Finally, the very sharp peak seen in Region 3 
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Fig. 13. Root-mean-square (rms) of the differences between the normalized flux of six model spectra (P, A1/A2, A3, Ml, M2/M3, 
M4) and the mean of these six spectra, averaged over 2 A-wide wavelength intervals. The different panels show four different 
wavelength regions, referred to in the text. 



in Fig. 13 close to 8810 A is due to deviations for a strong Mg I 
line in two of the models. The line is absent in Ml and narrower 
than in the others in P (see Fig. 14 i). A comparison with the 



observations supports the models with larger broadening. 

In general, the outcome of the comparison can be summa- 
rized in the following way: 

- Good agreement is seen for strong atomic lines, with a few 
exceptions (e.g. the wings of Mg I 8807 A). 

- Medium-weak atomic lines are in partial agreement. 

- Large differences are seen for the IR Ca II triplet lines. 

- Molecular lines differ to varying degrees. 

5. Discussion 

As pointed out above, there seems to be no obvious difference 
between the results from different 'code-families' (MARCS, 
ATLAS). Experiment 2, the 'fixed parameter comparison', 
revealed quite significant differences in the synthetic line 
strengths, in particular for some medium to weak lines, molecu- 
lar lines, and for the Call IR triplet. The main reason for these 
differences seems to be the different line parameters used by the 
different groups. However, the spectrum modelling differences 
are only one aspect leading to the resulting scatter of the derived 



parameters of, e.g.. Star 1. In the following we briefly discuss 
various aspects that may affect the result. 



5. 1 . Method of parameter determination 

In an ideal world, observations with their uncertainties would 
be combined with models and their uncertainties via a Bayesian 
approach to estimate stellar parameters, including correlations 
between parameters. This is especially difficult in the real world, 
because uncertainties are often difficult to estimate. Some of the 
observing parameters (e.g., spectral point-spread function, tel- 
luric correction) require complex models to determine uncer- 
tainties, and, to our knowledge, no modelling code currently 
propagates uncertainties in, e.g., oscillator strengths and pres- 
sure broadening coefficients through to wavelength-dependent 
uncertainties on the output spectra. Therefore, different groups 
tried different approaches to fitting what they regard as the most 
reliable parts of the spectra with differing weighting schemes. 

The participating teams have chosen a remarkable variety of 
approaches for deriving the stellar parameters. The groups can 
be divided into those which measured equivalent widths from the 
provided spectra which they compared to calculated equivalent 
widths (A2, A5, T), and those which employed spectrum synthe- 
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Fig. 14. Normalized flux differences from the mean for six different models, for four different wavelength regions. The six model 
codes appear across the bottom of panel (a) in the colours used in the plots, (a) blue region with relatively low rms difference; (b) 
region containing Ha; (c) region containing the IR Call triplet 8498 A line; (d) region with high rms difference due to differences 
in a strong Mg I line at 8807 A. 



sis. The latter groups synthesized either continuous wavelength 
intervals (M2, M3a, M4, M5, A3, P, C), or intervals around se- 
lected spectral lines (Ml, M3b, M6, Al, A4; see Fig. [T). The 
main difference between the two M3 analyses was in the char- 
acteristics of the spectra that were fed through an automated 
pipeline into the parameter determination code MATISSE. M3a 
was tailored to low resolution (R ~6500) spectra covering al- 
most continuous wavelengths within a region of ~300A. M3b 
was applied to spectra at a higher resolution (R ~15 000), using 
discontinuous wavelength regions totaling ~435A spread over 
~1800A. 

The M2, M5, M6, Al, A3, and T, analysis started by deter- 
mining approximate stellar parameters (or at least one param- 
eter) from a comparison of observed photometry and synthetic 
colours (see Tables [T2| to [TS] and Section [3.2.3| l. Ml used the 
colour information at a later step to partly lift the degeneracy be- 
tween Jeff/ log g and metallicity. The remaining groups did not 
use the provided colour information. 

The experiment participants used different methods to cal- 
culate equivalent widths or synthetic spectra. A2 applied auto- 
matic measurement of equivalent widths from synthetic spectra 
calculated with SYNTHE, using fixed- width integration around 
the center of each spectral line. A5 used the plane-parallel code 
MOOG (original version described in |Sneden|ri973p to calcu- 



late the equivalent widths, while T employed the method de- 
scribed in Tsuji ( 1986 1. Al, A3, and A4 worked with the plane- 
parallel line-formation code SYNTHE (Kurucz 1993b| l for spec- 
trum synthesis. Ml used a spherically-symmetric spectrum syn- 
thesis code developed for the study of line profile variability 
in AGB stars, described in |Nowotny et al. ( |2010 their Section 
2.4). M2, M3, and M5 used the spectral line formation code 
TURBOSPECTRUM ( [Alvarez & Plez|1998| l, while M4 used the 
code BSYN developed in Uppsala. Both TURBOSPECTRUM 
and BSYN provide the option to compute the radiative trans- 
fer in spherically-symmetric geometry, and both have a number 
of subroutines in common with the MARCS code. M6 used the 
MOOG code in spectrum-synthesis mode. The P and C analyses 
employed the spectrum synthesis integrated in the PHOENIX 
and CODEX packages, respectively. 

The methods used for parameter optimization with the pro- 
vided spectra are briefly summarized as follows. Ml, M2, M4, 
M5, A2, A3, P, and C based their fit on minimizing a weighted 
squared difference (x^) between the observed and the synthetic 
spectrum, either for the whole wavelength range available or for 
selected regions containing parameter-sensitive features. Some 
of the groups complemented the automatic fit with an inspection 
by eye (Ml, M2, M5). M3 used the MATISSE algorithm, which 
has been developed for fast automatic parameter determination 
of low- to high-resolution stellar spectra, and is based on a lo- 
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cal multi-linear regression method (for details see Section [3. 1.3| l. 
Al verified by eye that the observed spectrum was reproduced 
by the model spectrum with the photometric T^ff values, and ad- 
justed the log g and [Fe/H] values manually in order to fit se- 
lected features sensitive to these parameters. M6, A4, and A5 de- 
rived the stellar parameters by removing line abundance trends. 
This "classical" method effectively combines the information 
from spectrum points within spectral lines, and uses knowledge 
on spectral line behaviour to infer the parameters. The differ- 
ence was that M6 and A4 used spectrum fits for each line, and 
A4 also fitted the wings of strong lines, while A5 worked with 
equivalent widths. In the T analysis, T^ff was derived from pho- 
tometry, a fixed value was chosen for log g, and the abundance 
analysis was performed as described in Tsuji (2008 1. 



A prime issue that has to be discussed is how these different 
strategies affect the result. There are good reasons to start with 
photometry data to constrain the stellar parameters. Photometric 
colours have the advantage of providing information over a wide 
baseline in the spectrum. Some of them - of course depending 
on the spectral types of the target stars - are nearly independent 
of metallicity over a wide range of temperatures (for cool stars 
V -K can be used very well), while others involve a dependence 
on surface gravity. A clear disadvantage, however, that can intro- 
duce considerable error for this approach, is interstellar redden- 
ing, although this will not play a significant role in the context of 
our experiment. Narrow-band photometry - e.g. using filters to 
study the strength of individual molecular bands (see e.g. [Wing 
2007] l - may be a way to circumvent this weakness. 

Furthermore, photometric colours are affected by the actual 
chemical composition of the star, especially for cool objects with 
altered C/O ratios (e.g. Star 4 in Experiment 1). In such cases, 
the uncertainty in the derived Teff is expected to increase sig- 
nificantly. To make use of the photometry, a calibration relating 
colours and parameters is required, which can be either purely 
theoretical or semi-empirical. A considerable number of such 
calibrations have been published throughout the years. Thus, it 
is not surprising that each of the groups using the colour in- 
formation in Experiment 1 applied a different calibration (see 
Tables [12] to [TS] ). The resulting estimates for T^s agree rather 
well, except for Star 4, as expected from the discussion above 
(see also discussion in Section 3.1.5 1. For logg' and [Fe/H], the 



values disagree in the few applicable cases. 

All of the methodologies used here for the parameter de- 
termination suffer from several sources of uncertainty. First, in 
the squared-difference minimization method, similar values 
can be obtained for a certain range of models, and so the "best 
parameters" must be constrained by additional considerations. 
An interesting aspect is that the groups using this method de- 
fined the difference measure (x^) in different ways, in particu- 
lar regarding the weight factor It is not clear how this might 
affect the outcome of the analysis. Note that we did not com- 
pare individual values to assess the results from the differ- 
ent groups, as those will also depend on the length and location 
of the wavelength sections. Second, none of the methods pro- 
vides a straight-forward way for estimating the uncertainties in 
the derived parameters, except for MATISSE (M3) which auto- 
matically provides error estimates based on the signal-to-noise 
ratio. Third, separating the effects of different line-broadening 
parameters such as micro- and macro-turbulence or rotation may 
be difficult in the case of the spectral synthesis method. On the 
other hand, the full spectrum fit can - in principle - better handle 
spectra with lots of blended lines, and it includes a large number 
of spectral indicators at once. 



Since M2 and M3 used identical model spectra for their anal- 
ysis, the differences in the results are entirely due to the dif- 
ferent methods applied. In addition, the influence of different 
methods can be seen by comparing the M3a and M3b results. 
M2 and M3b both used the optical wavelength region (4900 - 
6800 A) for the fit, although slightly different subintervals were 
selected. Furthermore, M3b worked at decreased spectral resolu- 
tion, while M2 used the full resolution. The resulting Teff values 
are comparable for both a Tau and a Cet, while log g and [Fe/H] 
differ. This is in line with the findings of M4 that the optical re- 
gion is mostly sensitive to Teff variations. Hence, the "best-fit" 
values of log g and [Fe/H] will be rather arbitrary. Using only 
the RVS wavelength region (8470-8740 A) and a low resolu- 
tion, M3a gives different results for all three parameters. The Teff 
values are higher for both benchmark stars (compared to M3b), 
while log g and [Fe/H] are higher for one star and lower for the 
other. 

The P analysis for a Tau resulted in the highest value for 
loggf. This could be explained by low values of the pressure 
broadening parameters for strong atomic lines used by the P 
code, apparent in Fig. 14 c and d. To achieve a fit to broad ob- 
served line wings, the P analysis requires a higher log g value 
than the other models. 



5.2. Limitations in accessible parameter range 

Typically, models are optimized for application to a specific 
range of stellar parameters. In Experiment 1 we also see this 
in some limits of the tested parameter range. For example the 
CODEX code (C) has an upper limit in log g of 0.0, and thus the 
'best fit' result in the case of Star 2 is naturally somewhat offset 
from the findings of the other groups. 

most participating groups ex- 



As mentioned in Section 4.2 



plored only a small volume in parameter space around values 
estimated from photometry or comparisons with spectra of stan- 
dard stars. Only M3 and A2 based the analysis of Star 1 and 
Star 2 on a large model grid, and they obtained results similar to 
the others, except for T^^ in the case of A2. 

5.3. LIE and non-LTE 

As one of the underlying assumptions in this comparative work, 
all the participants used LTE stellar model atmospheres with 
LTE line formation (except for the M6 Ca abundance analysis 
described below). A full non-LTE treatment of all elements in 
the atmospheres of cool giants is a very complex task that can 
not yet be achieved due to the lack of quantitative atomic data. 
Nevertheless, non-LTE theoretical models of giant stars have 
been computed by Short & Hauschildt ( 2005 1, taking into ac- 
count the main electron contributors. They showed that non-LTE 
effects on the structure of the atmosphere increase with decreas- 
ing metallicity, especially deep and high up in the photosphere. 
The change in the structure will change the line formation do- 
mains and the intensity of the photospheric lines. Fortunately, 
the change in the structure seems to be negligible at solar metal- 
licity. Thus, one can use LTE model atmospheres for a partial 
non-LTE analysis by applying line synthesis of a given element 



in non-LTE. This has been done by many authors, e.g., Thevenin 
|& Idiart|([T999]l;|Shchukina & Trujillo Bueno| ( |200T) . More re- 
cently, Mashonkina et al. (201 1 1 showed with a very complete 
model atom of Fe I/II that the non-LTE corrections do not ex- 
ceed 0.1 dex for solar and mildly low metallicities. 
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The LTE assumption may strongly affect the surface gravity 
determination when this astrophysical parameter is constrained 
by the ionization equilibrium to follow Saha's law. The surface 
gravity deduced from the photometry and evolutionary tracks is 
usually used as a starting point for spectroscopic iteration be- 
tween two ionization stages of iron. This latter procedure relies 
on the assumption of LTE and tends to give LTE values of Teff 
and log g smaller than those from photometric determinations. 
Therefore, this can also impact the iron abundance analysis. A 
temporary solution to this problem can be to trust the surface 
gravity given by the photometric estimation, and to use only Fe II 
lines, which are known to be free of non-LTE effects (see, e.g. 
Mashonkina et al.t201 1 ), to determine [Fe/H]. An investigation 
of non-LTE effects on Fe line formation is not within the scope of 
the present paper However, the M6 team determined the abun- 
dance of Ca in a Tau using non-LTE radiative transfer in an LTE 
model atmosphere, as described in the following paragraph. 

Ca I/Call analysis of a Tau: The M6 team aimed to find the 
same abundance for Ca I and Ca 11. It was decided to compare the 
abundance determined using the IR triplet (8498, 8542, 8662) A 
with the abundance determined using the optical triplet (6102, 
6122, 6162) A. A non-LTE analysis was performed, using a Cal 
atom with 153 energy levels (taking into account fine structure), 
2120 bound-bound radiative transitions and 1 1476 bound-bound 
collisional transitions with electrons, and 8 1 photoionization ta- 
bles from TopBase. They used a Call model atom with 74 en- 
ergy levels, 422 bound-bound radiative transitions, 2628 bound- 
bound collisional transitions with electrons, and 40 photoioniza- 
tion tables from TopBase. The contribution of inelastic collisions 
with neutral hydrogen was neglected due to the lack of accurate 
quantum-mechanical calculations to determine their importance, 
and the Stark effect in the line-broadening parameters was also 
neglected (see [Merle et al.|201l| l. Spherical MARCS model at- 
mospheres with standard composition were adopted. The code 
MULTI2.2, slightly modified for the collisional transition part, 
solves the radiative transfer and statistical equilibrium equations 
consistently and computes line profiles, equivalent widths and 
contribution functions for all the radiative transitions. They per- 
formed several computations for the model atmosphere 4000/1/0 
(Teif/log g'/[Fe/H]) with different calcium abundances for Cal 
and Ca II and decided by eye the best metallicity fitting both the 
Cal lines and the Call triplet lines. This turns out to be [Ca/Fe] 
= -0.3 ±0.1 with an adopted [Fe/H]=0.0 (see Table[T2]i. We note 
that the core of the Ca II IR lines is still poorly reproduced, be- 
cause of the existence of a chromosphere where the line cores 
form (see |Berio et al.|201 l| l. 



5.4. Geometry 

In addition to Teff , log g and abundances, stellar spectra are in 
principle also affected by the extension of the atmosphere, which 
can be parameterized by (reff/5770)(/?/RQ)(M/Mor'. If the at- 
mospheric extension could be modelled from spectra, then this 
would uniquely constrain M, L and distance independently of 
evolutionary models. None of the models submitted in the study 
of stars 1 through 4 came up with a meaningful constraint on 
extension, but this by itself does not mean it is insignificant. 

Especially extended atmospheres can show very strong ef- 
fects of extension such as molecular shells, where it is not even 
obvious if a strong line in a cool extended shell will appear 

I. Extreme effects 
t to model com- 



pletely in studies such as this paper, because they necessarily re- 
quire departures from hydrostatic equilibrium (e.g. pulsation and 
shocks) in order to elevate the material. The effects in stars with 
only low-amplitude oscillations such as Star 2 are much smaller, 
but are significant enough that many modelling teams chose to 
use spherical codes. This mirrors the experience of the various 
groups that even in a purely hydrostatic world the sphericity ef- 
fects on spectra increase towards cooler giants. 

The extension of an atmospheric model is typically param- 
eterised by a single statement: whether the model is spherical 
(with a reasonable value of extension) or plane-parallel. For 
Teff greater than 4000 K, Heiter & Eriksson (2006) show typi- 



in absorption or emission, Ohnaka (e.g. 2004 



observed in highly evolved stars are difficu 



cal abundance differences between spherical and plane-parallel 
codes of 0.1 dex. This is comparable to or smaller than the metal- 
licity errors estimated from models in this paper, but suggests 
that if sphericity is not taken into account, the atmospheric pa- 
rameters will be mis-estimated. We suggest that further work 
on this topic is required, informed by detailed observations of a 
range of stars with well-known parameters and a range of atmo- 
spheric extensions. 

5.5. Atomic and molecular line data 

Accuracy and completeness of atomic and molecular data play 
a key role in fitting observed spectra with synthetic ones. In the 
experiments presented in this article we find various approaches 
to the compilation of input data. All MARCS-based analyses re- 
lied on atomic data extracted from the VALD database (in 2008- 
2009), while all ATLAS-based analyses used atomic data from 
the Kurucz web site. A few groups calibrated the gf/-values on 
high-resolution spectra of standard stars - M3a and M5 on the 
Sun and Arcturus, and A4 on several standard stars, including 
the Sun, Arcturus, and yu Leo. A2 verified the atomic line data 
for the cases of the Sun and Procyon. Molecular data were in- 
cluded in all calculations except for M6 and A5. Typically, line 
lists for 8-1 1 different species were included. In the case of M2, 
M3, and M4, the data stem from the same sources (described 
in 'Gustafsson et al. '2008). For Ml, the data sources are given 
in Lederer & Aringer (2009), and for M5 they are specified in 
Section |3.1.5| Al, A2, A3, A4, P, and C obtained molecular line 
data from the same sources (via the Kurucz website), see respec- 
tive subsections of Section 3 for details. Data sources for T are 
given in Tsuji (2008). 

The VALD database for atomic lines was therefore the prime 
source of atomic data, although the line lists have been extracted 
from different versions of the database. Sources for molecu- 
lar line lists are less homogeneous, and also the number of 
molecules included varies from group to group. All line lists, 
both for atoms and for molecules, still require significant im- 
provement. A way around missing line data is offered by the use 
of astrophysical oscillator strengths. 

However, the derivation of astrophysical gr/-values requires a 
reference object, for which atmospheric parameters and element 
abundances are known precisely. Currently, this is only the case 
for the Sun (even in that case some uncertainties remain). As the 
nature and strength of spectral lines appearing in the object under 
study (e.g. a cool giant) and the Sun will differ significantly, the 
astrophysical ^^/-values will be in most cases of limited use. 

A possible approach would be to work differentially using a 
"chain" of reference stars, starting from parameters close to so- 
lar and leading towards the required area in parameter space. We 
assume that for each neighbouring pair of stars it will be possible 
to define a common set of spectral lines. The line data calibrated 
on the star "closer" to the Sun will then define the element abun- 
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dances of the "new" star. However, astrophysical gi/-values will 
always remain an unsatisfying solution. Good experimental data 
are certainly highly desired. 

We note that oscillator strengths are only one element of 
the line data set required for each transition. Various damping 
constants used to parametrize coUisional processes are equally 
important and equally uncertain. Differential work can remedy 
these uncertainties to a certain extent, in a similar way as for 
gf/-values. For certain transitions, specialised theories are avail- 
able, as for example for the hydrogen line^^ but these are not 
implemented in all codes. 

5.6. Convection 

An important physical process to be included in the calculation 
of model stellar atmospheres is convection. The main effects of 
convection are 1) the contribution of the convective flux to the 
energy transport, leading to a modified temperature distribution 
compared to the purely radiative case; and 2) the generation of 
horizontal and vertical velocity fields, which modify the posi- 
tions and shapes of spectral lines. The efficiency of convective 
energy transport and its effect on the temperature structure in the 
line-forming region depends on all three atmospheric parame- 
ters (Teff, \ogg, [M/H]) in a complex way (the contribution of 
the convective flux to the total flux is shown e.g. in Fig. 4 in 
Heiter et al.|2002 l. At refF*4000 K, the efficiency of convection 
reaches a maximum and seems to be independent of gravity or 
metallicity. 

In all of the atmospheric models used for the experiments 
described in the present paper, the first effect of convection is in- 
cluded in an approximative way, based on the mixing length the- 
|1958) l. The practical implementation 



ory (MLT, Bohm-Vitense 



of the MLT varies somewhat between the different codes. In the 
ATLAS models, MLT convection is implemented as described 
in |Kurucz| ( |T9 70 ) and Castelli ( 1996), similar to the formulation 
according to Mihalas ( 1978 ). A variant of the latter is used in the 
PHOENIX models (Ludwi g et al.|[2002l). In the M ARCS mod- 



Henyey et al. ( |1965| l is used. This 



els, the MLT formulation of 

variation in the implementation results in different numbers of 
free parameters, and different meanings of these parameters, for 
each model, hampering a direct comparison (see Section 5.3.1 
in [Ludwig et al. 2002 for a discussion on the importance of the 
MLT formulation). 

The free parameter dominating any discussion of MLT con- 
vection is the "mixing length," i.e. the distance which convec- 
tive elements can travel before they dissolve. It is usually spec- 
ified in terms of the local pressure scale height, and typical val- 
ues are between 0.5 and 3, depending on stellar parameters and 
whether spectral features or global stellar parameters are mod- 
elled. However, a number of additional parameters have to be 
set. The size of convective elements is usually assumed to be 
a single value and equal to the mixing length. The need for 
this "one-eddy" approximation is overcome in the alternative 
"full-spectrum-turbulence" convection models (e.g. |Canuto &| 
Mazzitellil|1991| l, where the distribution of kinetic energy gen- 
erated by convection over different spatial scales is computed 
from a turbulence model. This treatment of convection has been 
implemented in ATLAS models by Heiter et al. ( 2002] l and was 
compared to MLT models for Teif >4000 K, log g>2 and a range 
of metallicities, but has not been used by any of the groups in the 
present paper. Additional free parameters are needed to specify 
the geometry of the convective elements and the temperature dis- 



tribution within them. This determines the energy lost by radia- 
tion and thus the efficiency of convection. Furthermore, viscous 
energy dissipation is parametrized by a coefficient in the equa- 
tion for the convective flux. The default values for these param- 
eters are different in the ATLAS and MARCS models. However, 
since all free parameters basically are scaling factors in a single 
equation, this may be compensated by using different values for 
the mixing length. 

An additional degree of freedom is provided in the MLT im- 
plementation of ATLAS models through the "overshooting op- 
tion." A smoothing procedure allows for a positive convective 
flux in stable layers next to the convection zone. Castelli et aLJ 
( |1997j investigated the effect of this option for reff>4000 K, 
loggf>2.5 and a range of metallicities. Comparison with obser- 
vation indicated that except for the Sun, the atmospheric param- 
eters derived with different methods were more consistent when 
the overshooting option was switched off. 

The second effect of convection, the three-dimensional ve- 
locity field, is not treated at all by MLT. Additional Doppler 
broadening of spectral lines caused by convective motions is 
taken into account in the one-dimensional models by specifying 
the ad-hoc parameters micro- and macroturbulence (in units of 
velocity, see modelling descriptions in Section |3] and Tables 12 
to[T5]). 

Realistic modelling of the effects of surface convection 
on stellar spectra can only be achieved with time-dependent, 
three-dimensional, radiation-hydrodynamic numerical simula- 
tions (hereafter referred to as "3D models"). Such simulations 
have been done for the Sun for several decades, and have become 
increasingly accurate. This has been demonstrated by compari- 
son to a variety of observations, not the least the profiles of nu- 
merous individual resolved spectral lines. Without applying any 
micro- or macroturbulence parameters, the wavelength shifts and 
asymmetries of the solar line profiles are in general very well re- 
produced (see, e.g., the detailed review by |Nordlund et al. 2009 1. 

3D surface convection models and spectra for stars other 
than the Sun are not yet widely available, due to the large 
amount of computing power and human post-processing time 
required. However, individual regions in the parameter space 
have been investigated thoroughly. The most relevant study in 
the context of this paper is that by |Collet et al.| ( |2007| , who 
examined the effects of convection on spectral lines for red gi- 
ant stars with log g - 2.2 at two different effective temperatures 
and [Fe/H] =0,-1,-2,-3. The largest effects when compared 
to ID model atmospheres where found at the lowest metallici- 
ties, while the effects were minor at solar metallicity (abundance 
differences of up to 0.8 dex and within ±0.1 dex, respectively). 
However, these models were significantly hotter and at lower 
log g than the stars investigated in the present paper 

The situation may improve rapidly in the near future, as at 
least two groups are computing larger grids of 3D models, us- 
ing two different simulation codes. The "CIFIST grid" ( [Ludwig | 
et al.|2009| l and the "StaggerGrid" ( [CoUet et al.|20n] ) both will 



http://www.astro.uu.se/<tilde>barklem/hlinop.html 



include on the order of 100 models for main-sequence and red 
giant stars (25-30 Teff/logg' combinations and four metallici- 
ties). Within the latter project, average temperature structures 
calculated from the full 3D models will be made publicly avail- 
able, as well as corresponding values for micro- and macrotur- 
bulence for implementation with classical line-formation codes. 
This will be the first step towards a more realistic spectrum anal- 
ysis of a variety of stars. We note that the parameters of the two 
benchmark stars a Tau and a Cet lie just outside of the parame- 
ter ranges investigated for convection effects so far, both in grids 
of ID models and in sets of 3D models. However, both of the 
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future 3D grids include a solar-metallicity model at (reff, log g) 
- (4000 K, 1.5), which should be close enough for establishing 
the importance of convection for spectrum modelling for these 
stars. 

Finally, we point out that the spatial extension of hydrody- 
namic simulations can be set up in two different ways. In the sur- 
face gravity domain of red giants, the "box-in-a-star" approach 
is applicable (cf the two model grids mentioned above). In this 
approach, the simulation is performed for a box inside the stel- 
lar atmosphere, covering several convective cells. The horizontal 
boundaries of the box are periodic, the vertical boundaries are 
open, and the extensions of the box are small compared to the 
stellar radius. However, the characteristic size of the convective 
cells scales with pressure scale height and thus increases with 
decreasing gravity. In the extreme case of red supergiants such 
as Betelgeuse, the atmosphere contains only a few giant convec- 
tive cells, and hydrodynamic simulations of the "star-in-a-box" 
type must be employed ( Freytag et al.||20(j2 i. a Tau and a Cet 
lie in the transition region between these two approaches (close 
to \ogg- 1, Chiavassa et al.|20TT i. Since the influence of con- 
vection on spectrum modelling has not been investigated in this 
region so far, these two stars are very interesting test objects for 
future studies (see, however. Section 3.5 in |Kucinskas et al.|2005| 
for a 3D hydrodynamical model and convection-related effects 
on colours of a late-type giant). 



5.7. Other issues 

An often neglected fact is that model atmosphere and synthetic 
spectrum codes rely on a wealth of physical input data besides 



the spectral line data discussed in Section 5.5 Apart from pa- 
rameters describing convective processes, an example for these 
are data used for chemical equilibrium calculations (establishing 
the equation of state), namely partition functions and molecular 
dissociation energies. Another important ingredient are data de- 
scribing bound-free and free-free transitions providing the con- 
tinuous opacity. Although widely used standard references ex- 
ist for much of these data (e.g. the "Kurucz routines" for con- 
tinuous opacities, or Irwin |[T98T for partition functions), these 



are evolving slowly and often are not described in detail. These 
"hidden data" leave room for differences of unknown degree be- 
tween codes from different authors. Furthermore, they can lead 
to inconsistencies between model atmospheres and spectral syn- 
thesis. In this regard, the first step towards consistent spectrum 
modelling is to use spectrum synthesis and model atmosphere 
codes developed by the same author or group. 

For example, the MARCS model atmosphere code and the 
BSYN synthetic spectrum code share many subroutines (e.g. for 
the chemical equilibrium) and the computation of the continu- 
ous opacities is consistent in both codes. The same is true for 
MARCS and TURBOSPECTRUM, and the "Kurucz suite" con- 
sisting of ATLAS9, ATLAS 12, and SYNTHE. 

6. Conclusions 

Estimating stellar parameters by comparison of observed and 
synthetic spectra is always affected by inaccuracies in line data 
and assumptions of model stellar atmospheres for cool stars (cf. 
the extensive discussion of this issue in the context of cool, 
evolved giants by Gustafsson 2QQ7\ . Due to these limitations 
even a 'perfect' fit of a spectrum cannot be seen as proof for 
a perfect model. We stress that it is not the intention of this pa- 
per to rank the various existing models of stellar atmospheres, 
but to learn from the effects caused by their differences. 



Experiment 1 clearly illustrates the need to be cautious when 
comparing or combining stellar parameters that were derived us- 
ing different model atmospheres and fitting strategies. We tried 
to illuminate various possible influences, but a clear trend in 
terms of a systematically higher or lower value for some pa- 
rameter for a given assumption could not be derived. This is a 
very complex problem with various aspects that seem to partly 
compensate for each other 

It would be desirable to repeat the experiment using differ- 
ent modelling codes, but exactly the same method (e.g. a chi- 
square fit to preselected wavelength regions). Also, differences 
in atomic and molecular line lists should be sorted out. It seems 
to be extremely difficult to implement the use of equal input line 
data in the various codes, simply because of different require- 
ments for the data format. This annoying obstacle could possibly 
be mitigated in the future by the use of a common infrastructure 
for atomic and molecular line data. An attempt to create such an 
infrastructure is currently underway in the VAMDC EU projecp] 
(Duberne t et al.|[2010| l. Finally, repeating the experiment using 
exactly the same tools may give a clearer idea about the reasons 
for different results between different codes. 

One could ask whether the data given for Experiment 1 are 
sufficient to derive a unique set of parameters. No doubt addi- 
tional data would set further constraints and resolve one or an- 
other degeneracy. A more extensive analysis of the existing data 
may produce some progress as well, but in our view it will be 
always limited by the accuracy of the line data and inaccuracies 
in our models, which definitely require further attention in the 
future. 

However, we think that Experiment 1 represents quite a typ- 
ical situation in observational astrophysics: a piece of the spec- 
trum and some photometric data are given, and some fixed strat- 
egy to analyze the data is adopted. With our limited experiment 
we want to point out that the atmospheric model and the anal- 
ysis method used can have a significant effect on the absolute 
values of the derived stellar parameters. Bringing together such 
a variety of groups fitting stellar spectra by the use of various 
methods, this paper provides also an unprecedented snapshot of 
the current status of this field. 

For future analysis projects, we recommend working differ- 
entially as far as possible. For stellar samples with parameters 
significantly different from the Sun, a set of suitable benchmark 
stars may serve as stepping stones for a piecewise differential 
analysis. 
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